FINAL  REPORT 


Improved  Understanding  of  In  Situ  Chemical  Oxidation 
Contaminant  Oxidation  Kinetics 


SERDP  Project  ER-1289 


December  2007 


Paul  Tratnyek 

Oregon  Health  &  Science  University 

Jamie  Powell 

Oregon  Health  &  Science  University 

Rachel  Waldemer 

Oregon  Health  &  Science  University 


Distribution  Statement  A:  Approved  for  Public  Release, 
Distribution  is  Unlimited 


t1  - SERDP 

Tdoe 

Strategic  Environmental  Research  and 
Dcve  lopmenl  P  ro  g  ra  m 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

DEC  2007  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2007  to  00-00-2007 

4.  TITLE  AND  SUBTITLE 

Improved  Understanding  of  In  Situ  Chemical  Oxidation  Contaminant 
Oxidation  Kinetics 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Oregon  Health  &  Science  University, Department  of  Environmental  and 
Biomolecular  Systems, 20000  NW  Walker  Rd,  Beaverton,  OR, 97006-8921 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

58 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


This  report  was  prepared  under  contract  to  the  Department  of  Defense  Strategic 
Environmental  Research  and  Development  Program  (SERDP).  The  publication  of  this 
report  does  not  indicate  endorsement  by  the  Department  of  Defense,  nor  should  the 
contents  be  construed  as  reflecting  the  official  policy  or  position  of  the  Department  of 
Defense.  Reference  herein  to  any  specific  commercial  product,  process,  or  service  by 
trade  name,  trademark,  manufacturer,  or  otherwise,  does  not  necessarily  constitute  or 
imply  its  endorsement,  recommendation,  or  favoring  by  the  Department  of  Defense. 


IMPROVED  UNDERSTANDING  OF  IN  SITU 
CHEMICAL  OXIDATION 

Technical  Objective  I:  Contaminant  Oxidation  Kinetics 


SERDP  Project  Number  ER-1289 


Principal  Investigator:  Paul  G.  Tratnyek 
Graduate  Students  and  Staff:  Rachel  Waldemer  and  Jamie  Powell 

Department  of  Environmental  and  Biomolecular  Systems 
Oregon  Health  &  Science  University 
20000  NW  Walker  Rd.,  Beaverton,  OR  97006-8921 

*Email:  tratnyek@ebs.ogi.edu 
Web:  http://cgr.ebs.ogi.edu 
Phone:  503-748-1023,  Fax:  503-748-1273 


OREGON 

HEALTH 


&SCIENCE 

UNIVERSITY 


FINAL  REPORT— OBJECTIVE  I 
May  2009 
Revision  1 


TABLE  OF  CONTENTS 


Page 

1 .  Acknowledgements . 1 

2.  Executive  Summary . 1 

3.  Project  Objectives . 3 

4.  Overview . 3 

5.  Outline  of  Tasks . 7 

6.  Determining  k"  for  Permanganate” . 9 

7.  Determining  K"  for  Activated  Hydrogen  Peroxide . 15 

8.  Determining  k"  for  Activated  Persulfate’ . 18 

9.  Correlation  Analysis . 25 

10.  IscoKin  database . 29 

11.  Oxidation  of  Chlorinated  Ethenes  by  Heat- Activated  Persulfate . 30 

12.  Effect  of  groundwater  constituents  on  ISCO  kinetics’ . 41 

13.  References . 45 


1 


LIST  of  TABLES 


Table  1: 
Table  2: 
Table  3: 
Table  4: 

Figure  1: 
Figure  2: 

Figure  3: 

Figure  4: 

Figure  5: 
Figure  1: 
Figure  2: 
Figure  3: 
Figure  4: 


Initial  conditions  in  the  Fenton’s  System  with  TCE 

k0bs  and  corresponding  activation  parameters  for  all  4  chlorinated  ethenes 

Arrhenius  parameters  for  chlorinated  ethenes 

£0bs  and  corresponding  activation  parameters  for  PCE  with  different 
concentrations  of  bicarbonate 


LIST  of  FIGURES 

General  scheme  summarizing  how  to  represent  the  disappearance  kinetics  of  a 
contaminant  of  concern  (COC)  due  to  reaction  with  an  oxidant  (Ox). 

Data  for  the  reactions  of  MnOf  with  3-chlorophenol  and  1,4-dioxane  (note  the 
difference  in  time  scales)  with  the  fit  to  a  simple  first-order  kinetic  model  (curve 
in  red)  and  Mn02-corrected  model  represented  by  equation  6  (curve  in  blue). 

Data  for  the  reactions  of  MnOf  with  3-chlorophenol  and  1,4-dioxane:  k0bs  vs. 
concentration  of  COC  plots;  the  slope  of  the  regression  line  is  used  to  determine 
k". 


Correlation  between  k"  obtained  from  pseudo-first-order  conditions  by  either 
analyzing  decreasing  concentrations  of  MnOf  in  the  presence  of  an  excess 
amount  of  COC  (this  study)  or  analyzing  decreasing  concentrations  of  COC  in  the 
presence  of  an  excess  amount  of  MnOf  ( k ”  data  obtained  from  the  literature). 

Summary  of  second-order  rate  constants  ( k ")  for  MnOf  and  compounds  in  all 
classes. 

Zero-headspace  reactor  used  for  to  measure  degradation  kinetics  by  iron-activated 
hydrogen  peroxide. 

Disappearance  kinetics  for  chlorinated  solvents  in  a  homogenous  Fenton  system: 
(A)  TCE,  and  (B)  1,1 -DCE. 

Modeling  of  kinetics  for  chlorinated  solvents  in  a  homogenous  Fenton  system: 

(A)  TCE,  and  (B)  1,1 -DCE. 

Typical  competition  kinetic  plots.  (A)  Change  in  absorbance  at  440  nm 
(indicating  change  in  concentration  of  PND  A  over  90  min);  (B)  Linear  version  of 
1st  six  data  points,  used  to  calculate  k". 


ii 


Figure  5: 
Figure  6: 

Figure  7: 
Figure  8: 
Figure  9: 
Figure  10: 

Figure  1 1 : 

Figure  12: 
Figure  13: 
Figure  14: 

Figure  15: 
Figure  16: 

Figure  17: 
Figure  18: 

Figure  19: 
Figure  20: 


LIST  of  FIGURES,  cont. 

Data  for  perchloroethene  (PCE)  does  not  conform  to  the  S-shaped  curve  typical  of 
competition  kinetics. 

The  values  of  k"  obtained  by  this  study  with  competition  kinetics  compared  with 
values  of  k"  for  'OH  and  SO4’  obtained  with  either  COC  or  SO4'  in  excess 
(literature  values). 

Summary  of  second-order  rate  constants  ( k ")  for  SO/-  and  aromatic  compounds. 

Summary  of  second-order  rate  constants  ( k ")  for  SO4'  and  aliphatic  compounds. 

Correlation  of  k"  for  oxidation  reactions  of  COCs  by  SCV”  to  E'homo- 

Scatter  plot  matrix  of  log  k"  of  all  compounds  obtained  from  experiments  or 
literature. 

Matrix  of  scatter  plots  of  log  k"  for  M  n(V  vs.  the  number  of  chlorine  atoms 
attached  to  the  double  bond  and  vs.  log  k"  for  O3. 

Screen  shot  of  the  initial  search  window  for  the  IscoKin  database. 

Pseudo-first-order  disappearance  of  PCE  at  five  temperatures. 

Concentration  over  extended  time  periods,  showing  non-pseudo  first-order 
behavior  at  longer  experimental  times. 

Arrhenius  plots  for  the  chlorinated  ethenes. 

Concentration  of  cis-DCE  and  trans-DCE  vs.  time  for  A)  the  reaction  of  cis-DCE 
with  persulfate  and  B)  the  reaction  of  trans-DCE  with  persulfate. 

Comparison  of  three  degradation  processes  for  PCE  as  a  function  of  temperature. 

Comparison  of  three  degradation  processes  as  a  function  of  temperature  for  TCE, 
cis-DCE,  and  trans-DCE. 

Effect  of  bicarbonate  on  the  oxidation  of  PCE  by  persulfate. 

Effect  of  bicarbonate  concentration  on  k0 bs  for  PCE  oxidation  by  persulfate 
activated  at  three  temperatures. 


iii 


LIST  of  ACRONYMS 


A 

BTEX 

°C 

CHCI3 

cr 

coc 

DCE 

DDT 

El/2 

Ea 

Eq\p 

Ehomo 

H202 

HMX 

HOMO 

IP 

ISCO 

ISCOKIN 

ISTR 

k0bs  or  k' 
k" 

LUMO 

M 

mg/L 

mM 

MnOf 

M11O2 

MTBE 

NAPL 

NOD 

•OH 

Ox 

O3 

PCE 

HSO5 

PNDA 

QSARs 

R 

RDX 

s2o82- 

S04" 

T 

TCE 

TNT 


pre-exponential  factor 
benzene,  toluene,  xylene 
degrees  Celcius 
chloroform 
chlorine  radical 
contaminant  of  concern 

1 . 1  -dichloroethene 

1.1.1  -trichloro-2,2-bis(p-chlorophenyl)ethane 
half-wave  potential 

apparent  activation  energy 

the  energy  difference  between  the  HOMO  and  the  LUMO 

energy  of  the  highest  occupied  molecular  orbital 

hydrogen  peroxide 

cyclotetramethylene-tetranitramine 

highest  occupied  molecular  orbital 

ionization  potential 

in  situ  chemical  oxidation 

database  containing  our  collection  of  k”  data 

in  situ  thermal  remediation 

pseudo-first-order  rate  constants 

second-order  rate  constants 

lowest  unoccupied  molecular  orbital 

molar 

milligrams  per  liter 
millimolar 
permanganate 
manganese  dioxide 
methyl  tert  butyl  ether 
non-aqueous  phase  liquid 
natural  oxidant  demand 
hydroxide  radical 
oxidant 
ozone 

perchloroethylene 

peroxymonosulfate 

p-nitrosodimethylaniline 

quantitative  structure-activity  relationships 

universal  gas  constant 

cyclotrimethylenetrinitramine 

Persulfate 

sulfate  radical 

absolute  temperature 

trichloroethylene 

trinitrotoluene 


1.  Acknowledgements 


The  majority  of  the  experimental  work,  and  much  of  the  data  mining  and  literature 
analysis  that  comprised  this  project  was  performed  by  Rachel  H.  Waldemer.  Additional 
contributions  were  made  by  Jamie  S.  Powell  and  James  T.  Nurmi.  Richard  L.  Johnson 
contributed  to  the  work  on  activated  persulfate.  We  thank  Eric  Hood  and  Neil  Thomson  for 
helpful  comments  on  some  of  documents  on  permanganate. 


2.  Executive  Summary 

The  overall  goal  of  Technical  Objective  1:  Contaminant  Oxidation  Kinetics  for  this 
project  was  to  provide  a  comprehensive  perspective  on  the  kinetics  of  oxidation  of  groundwater 
contaminants  by  ISCO  oxidants.  We  did  this  by  summarizing  previously  published  (second 
order)  rate  constants  for  all  of  the  putative  oxidants  (mainly  permanganate  anion,  hydroxyl 
radical,  and  sulfate  radical).  Where  important  data  gaps  are  discovered,  we  measured  rate 
constants  in  the  laboratory.  The  overall  dataset  was  used  to  develop  quantitative  structure- 
activity  relationships  (QSARs)  for  predicting  rate  constants  were  none  have  been  measured. 

First,  we  developed  a  novel  method  for  measuring  rate  constants  of  contaminant 
oxidation  by  permanganate.  Permanganate  disappearance  kinetics  are  pseudo-first-order  as  long 
as  contaminant  concentrations  are  in  excess.  The  colloidal  manganese  dioxide  that  forms  from 
permanganate  oxidation  gives  Beer’s  law  behavior,  and  a  procedure  was  developed  to  correct  for 
any  interference  this  caused. 

Using  this  method,  the  oxidation  of  most  contaminants  by  permanganate  was  shown  to  be 
first-order  in  contaminant  concentration  (second-order  overall)  by  varying  initial  concentrations 
of  contaminant.  The  second-order  rate  constants  obtained  with  our  method  are  in  good  agreement 
with  those  previously  reported  values  for  perchloroethylene  (PCE)  and  trichloroethylene  (TCE). 
Further  measurements  with  a  chlorinated  ethane  (1,1,1-trichloroethane)  and  a  BTEX  constituent 
(toluene)  show  that  reaction  of  these  contaminants  with  permanganate  is  slow.  Other  classes  of 
contaminants  that  were  tested  include  chlorinated  methanes  (carbon  tetrachloride,  etc.), 
explosives  (TNT,  RDX),  and  oxygenates  (MTBE). 

Activated  persulfate  produces  sulfate  radical,  hydroxyl  radicals,  and  other  reactive 
intermediates.  We  investigated  three  methods  for  activation  of  persulfate  (light,  heat,  and 
chelated  ferrous  iron)  and  showed  that  these  methods  produced  very  different  amounts  of  radical 
intermediates.  We  developed  a  method  for  measuring  rate  constants  for  contaminant  oxidation 
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by  sulfate  radical  and  applied  it  to  a  variety  of  contaminants  including  the  BTEX  compounds, 
chlorinated  methanes,  1,4-dioxane,  and  MTBE.  Rate  constants  for  reactions  with  sulfate  radical 
were  related  to  rate  constants  for  hydroxyl  radical  to  shed  some  light  on  the  relative  significance 
of  these  two  oxidants  during  ISCO  with  activated  persulfate  (under  laboratory  and  field 
conditions). 

While  chlorinated  ethenes  have  been  successfully  remediated  with  both  in  situ  chemical 
oxidation  (ISCO)  and  in  situ  thermal  remediation  (ISTR),  there  is  evidence  that  combining  these 
two  treatment  approaches  may  result  in  synergistic  advantages.  Sodium  persulfate  is  an  ideal 
oxidant  for  combination  with  ISTR  because  persulfate  exposure  to  high  temperatures  leads  to  the 
formation  of  highly  reactive  sulfate  radicals  in  addition  to  the  higher  reaction  rates  normally 
induced  by  higher  temperatures.  We  have  found  that  the  reactions  of  chlorinated  ethenes  and 
heat-activated  sodium  persulfate  are  first-order  with  respect  to  chlorinated  ethenes  for  at  least  the 
first  one  to  two  half-lives,  and  with  sodium  persulfate  in  excess,  the  resulting  pseudo-first-order 
reactions  can  be  described  by  the  Arrhenius  equation.  The  dichloroethenes  appear  to  mineralize 
more  completely  than  TCE  and  PCE  based  on  chlorine  recovery  analysis. 

Using  Arrhenius  parameters  (both  published  and  those  determined  here)  for 
permanganate  oxidation,  persulfate  oxidation,  and  hydrolysis  (the  ubiquitous  degradation  process 
that  could  be  significant  at  higher  temperatures),  comparisons  for  these  three  degradation 
processes  were  made  over  the  temperature  range  of  0  to  100  C.  It  was  determined  that  hydrolysis 
is  never  important  for  the  chlorinated  ethenes.  For  PCE,  permanganate  oxidation  may  be  more 
favorable  at  lower  temperatures,  but  persulfate  oxidation  is  more  favorable  at  higher  temperature. 

In  situ  chemical  oxidation  of  chlorinated  solvents  with  activated  persulfate  has  received 
considerable  study  in  the  last  few  years,  with  much  of  this  work  focused  on  comparison  of 
alternative  methods  for  activation  (heat,  iron,  peroxide,  etc.).  By  comparison,  data  on  the 
oxidation  of  other  contaminants  by  activated  persulfate  are  relatively  scarce.  We  revisited  the 
oxidation  of  explosives  with  heat-activated  persulfate  and  found  fairly  rapid  degradation  of  TNT, 
RDX,  HMX,  and  4-nitrophenol.  The  kinetics  generally  are  first-order  and  the  effect  of 
temperature  follows  the  Arrhenius  model.  The  reaction  mechanism  probably  is  not  dominated  by 
electron-transfer  to  sulfate  radical — as  is  the  case  with  most  organic  oxidations  by  activated 
persulfate — because  the  explosives  are  very  poor  electron  donors.  Instead  the  reaction 
presumably  involves  sulfate  radical  addition  or  substitution,  but  this  cannot  yet  be  confirmed 
because  no  products  have  been  identified  to  date. 
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3.  Project  Objectives 

The  use  of  in  situ  chemical  oxidation  (ISCO)  for  treatment  of  chlorinated  solvent  source 
areas  is  rapidly  increasing  as  Department  of  Defense  (DoD)  and  other  stakeholders  search  for 
remedial  approaches  that  reduce  long-term  operations  and  maintenance  requirements.  While 
ISCO  is  a  promising  technology  for  some  chemicals,  there  remains  significant  data  needs  related 
to:  i)  reaction  kinetics  for  common  DoD  contaminants;  ii)  the  effects  of  natural  oxidant  demand 
on  oxidant  mobility  and  delivery  under  varying  site  conditions;  and  iii)  the  effects  of  ISCO  on 
long-term  groundwater  quality.  This  study,  Strategic  Environmental  Research  and  Development 
Program  (SERDP)  research  project  CU-1289  (An  Improved  Understanding  of  In  Situ  Chemical 
Oxidation),  focused  on  addressing  these  data  needs  for  permanganate,  persulfate  and  Fenton’s 
reagent.  The  overall  goal  of  this  research  program  was  to  address  critical  research  needs  for  the 
improved  implementation  of  ISCO  using  Fenton’s  reagent,  persulfate  or  permanganate.  These 
critical  research  needs  include:  Technical  Objective  1,  the  development  of  a  comprehensive 
kinetic  perspective  on  the  kinetics  of  oxidation  of  common  DoD  contaminants  by  the  most 
commonly  used  oxidants,  permanganate  (MnOf),  Fenton's  reagent  (H202/Fe2+)  and  persulfate; 
and  Technical  Objective  2,  assess  how  soil  properties  (e.g.,  soil  mineralogy,  natural  carbon 
content)  affect  oxidant  mobility  and  stability  in  the  subsurface,  and  develop  a  standardized 
natural  oxidant  demand  (NOD)  measurement  protocol. 

This  project  report  is  divided  into  two  Final  Reports  that  address  each  of  these  critical 
research  needs:  Contaminant  Oxidation  Kinetics  and  Soil  Reactivity.  This  Final  Report  addresses 
Contaminant  Oxidation  Kinetics  and  was  completed  by  the  Oregon  Health  and  Sciences 
University.  The  Final  Report  for  Soil  Reactivity  was  completed  by  the  University  of  Waterloo. 
These  reports  build  upon  the  research  results  previously  presented  in  this  project’s  2003  and 
2004  Annual  Reports  submitted  to  SERDP  (Geosyntec,  OHSU,  and  UW,  2003  and  2005). 

4.  Overview 

In  the  last  few  years,  a  few  detailed  studies  have  been  reported  on  the  kinetics  and 
mechanism  by  which  permanganate  oxidizes  key  groundwater  contaminants,  mainly  the 
chlorinated  ethenes  (Huang  et  al.,  1999;  Hood  et  al.,  2000;  Yan  and  Schwartz,  2000).  The 
literature  on  oxidation  of  organics  by  Fenton’s  reagent  is  considerably  larger,  but  even  there,  the 
scope  of  most  detailed  studies  of  kinetics  and  mechanisms  tends  to  be  narrow  with  respect  to  the 
number  and  type  of  contaminants  considered  (e.g.,  Lipczynska-Kochany,  1991;  Sedlak  and 
Andren,  1991;  Leung  et  al.,  1992;  Huston  and  Pignatello,  1996;  Li  et  al.,  1997;  Tang  and  Tassos, 
1997;  Augusti  et  al.,  1998;  Bier  et  al.,  1999).  A  distinguishing  feature  of  this  project  was  that  it 
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focused  on  identifying  the  kinetic  and  mechanistic  characteristics  that  are  common  to  a  range  of 
contaminant  oxidation  reactions,  with  the  ultimate  goal  of  producing  a  comprehensive 
perspective  on  the  chemistry  of  ISCO  technologies. 

The  lack  of  a  comprehensive  perspective  on  the  chemistry  of  ISCO  technologies  is 
evident  from  the  numerous  inconsistencies  in  how  oxidation  by  permanganate  or  Fenton’s 
reagant  are  treated  by  most  environmental  chemists  and  engineers.  For  example,  the  strength  of 
both  oxidants  is  frequently  equated  with  their  very  positive  oxidation  potentials  by  advocates  of 
ISCO  technologies,  even  though  it  is  unusual  for  either  oxidant  to  react  by  electron  transfer 
(Tratnyek  and  Macalady,  2000).  For  permanganate,  the  only  detailed  mechanism  for  oxidation  of 
organics  that  has  been  described  in  the  environmental  literature  involves  electrophilic  addition 
and  rearrangement  of  chlorinated  ethenes  (Yan  and  Schwartz,  2000).  But  even  this  may  not  be 
the  whole  story,  as  recent  evidence  in  the  chemical  literature  has  shown  that  permanganate  reacts 
with  most  organics  by  hydride  abstraction  (Gardner  and  Mayer,  1995;  Mayer,  1998). 

With  respect  to  groundwater  applications  of  activated  hydrogen  peroxide  processes  (i.e., 
the  Fenton  reaction),  similar  inconsistencies  exist.  For  example,  it  is  almost  universally  assumed 
that  oxidation  is  mediated  by  hydroxyl  radical,  even  though  there  is  still  a  raging  debate  in  the 
chemical  literature  regarding  the  relative  importance  of  hydroxyl  radical  vs.  other  intermediates 
in  activated  hydrogen  peroxide  systems  (Sawyer  et  al.,  1996;  MacFaul  et  al.,  1998;  Walling, 

1998;  Goldstein  and  Meyerstein,  1999).  Remarkably,  the  reactive  intermediates  other  than 
hydroxyl  radical  that  have  received  the  most  attention  in  the  environmental  literature  are 
superoxide  and  hydroperoxide  anion  (Watts  et  al.,  1999),  and  these  species  are  hardly  mentioned 
in  the  debate  on  this  topic  among  the  chemists.  Clearly,  there  still  is  considerable  uncertainty 
over  how  activated  hydrogen  peroxide  degrades  organic  substance  in  water,  despite  many  years 
of  research  in  this  area. 

Despite  these  uncertainties  in  the  mechanism  of  contaminant  degradation  by 
permanganate  and  Fenton’s  reagent,  some  very  complex  and  detailed  kinetic  models  have  been 
described  for  these  oxidants  (e.g.,  Hong  et  al.,  1996;  Yan  and  Schwartz,  2000;  Zhang  and 
Schwartz,  2000).  The  development  and  application  of  these  models  has  advanced  our 
understanding  of  advanced  oxidation  processes  considerably,  but  they  have  rarely  been  applied 
directly  to  the  design  ISCO  applications  in  the  field  by  practitioners  of  these  technologies. 
Another  goal  of  this  project  is  to  bridge  the  gap  between  the  most  detailed  kinetics  models  and 
the  highly  simplified  kinetic  models  that  are  (often  implicitly)  used  in  engineering  design. 

The  philosophy  that  we  bring  to  this  problem  is  exemplified  in  the  influential  work  of  Prof. 
Jurg  Hoigne  and  coworkers  for  describing  the  kinetics  of  contaminant  oxidation  in  disinfection 
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and  photolysis  (e.g.,  Hoigne  and  Bader,  1983a;  Hoigne,  1990;  Tratnyek  and  Hoigne,  1991; 
Hoigne  and  Bader,  1993).  The  approach  is  basically  to  focus  on  parameterizing  the  model  as  the 
sum  of  independent  second-order  rate  laws  involving  measurable  reactant  concentrations. 
Although  there  is  nothing  novel  about  this  approach  per  se,  a  significant  challenge  lies  in 
codifying  this  process  at  a  level  that  is  rigorous  (in  a  physical  chemical  sense)  and  yet  still 
applicable  to  complex,  real-world  situations.  One  way  to  do  this  is  summarized  in  Figure  1. 

The  approach  taken  in  this  project  was  primarily  to  proceed  from  the  top  of  Figure  1  to  the 
bottom,  via  equations  1-4,  in  order  to  derive  generally-useful  properties  of  each  combination  of  a 
contaminant  of  concern  (COC)  and  oxidant  (Ox).  This  is  relatively  straight  forward  with 
permanganate  because  experiments  with  permanganate  are  formulated  with  the  oxidant  in 
excess,  such  that  [Ox]  is  known  and  constant.  In  contrast,  with  activated  hydrogen  peroxide  the 
primary  oxidant  is  hydroxyl  radical  (*OH),  which  is  a  reactive  intermediate  that  is  generated  at  a 
very  low  but  steady  state  concentration  (thus  the  subscript  “ss”).  A  similar  situation  exists  with 
activated  persulfate,  where  the  primary  oxidant  is  sulfate  radical  (SO/-).  Note  that  other 
experimental  variables  can  be  accommodated  in  this  scheme,  if  necessary.  Most  notably, 
temperature  can  be  described  by  inserting  the  Arrhenius  equation  between  equations  2  and  4. 
Once  the  fundamental  kinetic  data  (second  order  rate  constants,  k ")  are  available,  Figure  1  can  be 
applied  from  the  bottom  up  by  practitioners  to  estimate  the  performance  of  a  ISCO  treatment 
under  a  particular  set  of  conditions. 

In  light  of  the  importance  of  good  kinetic  data,  we  compiled  the  values  of  k”  from  the 
literature  for  oxidation  of  a  variety  of  contaminants  by  permanganate  (MnOf);  hydroxyl  radical 
(•OH),  often  thought  to  be  the  reactive  species  formed  by  the  Fenton  reaction  (Pignatello  et  al., 
2006);  and  sulfate  radical  (SO4 '),  thought  to  be  the  reactive  species  formed  by  activated 
persulfate  (Kolthoff  et  al.,  1951;  Peyton,  1993;  Couttenye  et  al.,  2002).  We  then  determined  the 
priority  data  gaps,  designed  several  methods  for  filling  the  data  gaps  for  these  oxidants,  and 
applied  these  methods  to  many  priority  COCs  (focusing  on  emerging  contaminants, 
contaminants  that  are  prevalent  at  military  sites,  and  contaminants  that  form  non-aqueous 
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Figure  21 .  General  scheme  summarizing  how  to  represent  the  disappearance  kinetics  of  a 
contaminant  of  concern  (COC)  due  to  reaction  with  an  oxidant  (Ox).  Pseudo  first  order  rate 
constants  ( k ')  are  used  to  describe  the  kinetics  of  a  particular  reaction  under  specific  conditions, 
from  which  second-order  rate  constants  (k")  are  derived  to  characterize  the  reaction  in  general. 
The  effect  of  pH  is  shown  where  it  involves  speciation  of  the  COC  into  protonated  and 
deprotonated  forms  (HA  and  A'),  but  an  analogous  treatment  is  possible  if  pH  effects  are  due  to 
speciation  of  the  Ox.  QSARs  that  can  predict  k"  can  be  expressed  as  linear  models  with  one  or 
several  descriptor  variables  (D)  and  the  associated  slopes  and  intercept  and  /%).  A  probe 
substance  can  be  used  to  determine  the  steady-state  (ss)  concentration  of  intermediate 
oxidants  such  as  hydroxyl  radical. 

phase  liquids  in  the  environment).  Then,  combining  new  and  previously  published  data  for  k", 
we  performed  correlation  analysis  to:  (/')  assess  the  overall  quality  of  the  data  set,  (ii)  investigate 
the  possibility  that  quantitative  structure-activity  relationships  (QSARs)  could  be  developed  to 
predict  values  of  k"  for  COCs  that  have  not  been  studied,  and  (Hi)  examine  whether  the  different 
oxidants  showed  similar  patterns  of  reactivity  with  the  same  set  of  contaminants,  indicating 
whether  different  oxidants  had  similar  mechanisms  of  reactivity. 

While  most  of  our  research  for  this  project  focused  on  compiling  and  generating  new 
values  of  k”,  there  were  some  situations — mainly  involving  activated  peroxide  and  persulfate — 
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for  which  it  was  more  practical  to  focus  on  pseudo-first-order  rate  constants  ( k0bs  or  k').  For 
example,  this  was  the  case  when  we  were  studying  the  effects  of  temperature  on  the  oxidation  of 
chlorinated  ethenes  by  heat-activated  persulfate  (Section  11).  Using  k0bs  as  our  kinetic  parameter 
allowed  us  to  measure  the  global  effect  of  temperature  on  the  many  reactions  that  are  generated 
once  persulfate  is  activated.  While  this  is  a  less  rigorous  approach,  studying  the  effect  of 
temperature  on  each  reaction  would  be  a  lengthy  process  and  was  not  possible  within  the  time- 
scale  of  the  study;  furthermore,  it  is  not  likely  that  the  detail  gained  from  this  approach  would 
change  the  conclusions  that  were  drawn  by  using  k0bs. 

5.  Outline  of  Tasks 

The  kinetics  part  of  this  project  was  structured  in  six  tasks:  (/')  assemble  and  summarize 
existing  kinetic  data,  (zz)  determine  and  prioritize  key  data  gaps,  (in)  validate  a  novel 
experimental  method  for  measuring  new  kinetic  data,  (zv)  apply  our  experimental  method  to  fill 
the  most  important  data  gaps,  (v)  perform  correlation  analysis  on  new  and  old  data  together  to 
validate  them  all  and  derive  models  that  can  predict  data  that  are  unavailable,  and  (vz)  begin 
calibration  of  a  quantitative  kinetic  that  can  describe  the  effect  of  competition  among  multiple 
oxidants.  Tasks  1-4  apply  separately  to  permanganate,  activated  peroxide,  and  activated 
persulfate,  whereas  Tasks  5-6  apply  to  all  the  ISCO  chemistries  collectively. 

Task  1.  A  exhaustive  literature  review  was  conducted  for  previously  published  kinetic 
data  on  the  oxidation  of  NAPL  constituents  by  permanganate,  hydroxyl  radical,  and  sulfate 
radical.  The  results  of  this  work  have  been  (or  will  be)  published  in  print  and  online,  and  are 
incorporated  in  Sections  6,  8,  9,  and  1 1  as  well  as  in  the  ISCOKIN  database  (Section  10). 

Task  2.  To  identify  and  prioritize  key  data  gaps,  we  summarized  the  data  from  Task  1 
into  a  variety  of  formats.  Figures  and  tables  of  these  results  were  presented  at  various  meetings 
and  the  first  annual  report  for  this  project.  Beyond  that,  they  have  not  been  published  because 
there  purpose  is  largely  help  with  prioritization  of  the  work  on  the  tasks  that  follow. 

Task  3.  In  preparation  for  experimental  work  to  fill  some  of  the  data  gaps  identified  in 
Task  2,  we  investigated,  evaluated,  and  validated  a  variety  of  experimental  approaches:  (z)  serial 
sampling  of  shaken  bottles  (Tratnyek,  1998;  Hood  et  al.,  2000;  Yan  and  Schwartz,  2000),  (zz) 
serial  sampling  and  continuous  monitoring  of  stirred  reactors  (Yan  and  Schwartz,  1999),  (z'z'z) 
serial  sampling  of  a  stirred  syringe  reactor  (Huang  et  ah,  1999),  and  (zv)  continuous  monitoring 
with  a  stopped  flow  reactor  (Nam  and  Tratnyek,  unpublished  data).  Different  methods  were  used 
for  different  oxidants,  so  they  are  described  below  together  with  the  results  for  each  oxidant. 
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Task  4.  Using  the  methods  developed  under  in  Task  3,  we  measured  new  values  of  k"  for 
permanganate  and  activated  persulfate  but  not  activated  peroxide  (because  the  existing  data  for 
•OH  reaction  with  organic  contaminants  is  very  extensive  and  extraction  of  k"  from  data  for 
Fe(II)  activated  peroxide  would  have  been  very  complex  and  time  consuming).  For 
permanganate  and  activated  persulfate,  we  determined  how  the  kinetics  of  contaminant  oxidation 
were  affected  by  key  experimental  variables  (e.g.,  temperature  and  pH).  These  results  are 
described  in  Sections  6-8,  11-12. 

Task  5.  The  combination  of  new  kinetic  data  from  Task  4  and  previously  published 
kinetic  data  from  Tasks  1-2  were  subjected  to  a  comprehensive  “correlation  analysis”  (Tratnyek, 
1998;  Farrell  and  Luo,  2002).  Correlation  analysis  is  a  very  powerful  tool  for  two  more  specific 
ends:  (i)  verification  that  the  data  are  free  of  erroneous  or  inconsistent  values,  and  (ii)  derivation 
of  quantitative  structure-activity  relationships  (QSARs)  for  predicting  data  that  are  not  available. 
The  results  of  this  are  summarized  in  Section  9;  they  have  not  yet  been  published  but  we  expect 
they  will  be  once  we  have  found  time  to  resolve  several  residual  issues  with  this  analysis. 

Task  6.  Although  the  kinetic  scheme  set  out  in  Figure  1  should  be  sufficient  to 
characterize  the  major  “intrinsic”  factors  that  determine  the  kinetics  of  oxidation  in  solution, 
there  are  a  host  of  environmental  factors  that  may  decrease  the  rate  of  contaminant  oxidation 
through  competition  for  the  oxidant.  These  competitors  for  oxidant  can  be  other  organic 
contaminants,  natural  organic  matter,  various  inorganic  solutes  like  carbonate,  and  mineral 
surfaces.  We  performed  experiments  to  characterize  the  competitive  effect  of  the  most  important 
of  these  groundwater  constituents,  and  the  results  are  included  below  in  Section  12. 


6.  Determining  k"  for  Permanganate1 ,2,3 

To  measure  k"  for  contaminant  oxidation  by  permanganate,  we  used  absorbance 
spectrometry,  an  efficient  method  that  allowed  us  to  measure  the  kinetics  of  oxidation  for  24 
contaminants — many  for  which  data  were  not  previously  available.  A  factor  that  complicates  the 
use  of  absorbance  spectrometry  to  follow  oxidation  reactions  involving  MnOT  is  that  the 
reduction  product,  colloidal  MnC>2,  contributes  to  the  absorbance  measured  at  525  nm  (an 
absorbance  maximum  for  M11O4")  (Mata-Perez  and  Perez-Benito,  1985;  Perez-Benito  and  Arias, 
1992;  Gardner,  1996;  Yan  and  Schwartz,  1999).  Fortunately,  the  measured  absorbance  of 
colloidal  MnC>2  suspensions  obeys  Beer’s  Law  until  the  particles  begin  to  flocculate,  resulting  in 
increased  scattering  of  light  by  the  larger  aggregates  (Mata-Perez  and  Perez-Benito,  1985;  Perez- 
Benito  and  Arias,  1992).  Our  methods  to  determine  where  Beer’s  Law  was  obeyed  in  our  system 
are  discussed  elsewhere  (Waldemer,  2004;  Waldemer  and  Tratnyek,  2004;  Waldemer  and 
Tratnyek,  2006). 

After  excluding  the  data  where  scattering  by  colloidal  Mn02  was  significant  (causing 
deviation  from  Beer’s  Law),  we  found  that  most  COCs  gave  data  for  absorbance  at  525  nm  vs. 
time  that  fit  pseudo-first-order  kinetics  with  a  modest  deviation  from  this  kinetic  model  as  the 
reaction  approached  completion  (Figure  2).  We  concluded  that  the  latter  effect  was  due  to 
increasing  light  absorption  by  suspended  MnC>2,  based  on  comparison  with  previously  published 
spectra  for  colloidal  MnC>2  suspensions  (Perez-Benito  and  Arias,  1992).  To  obtain  pseudo-first- 
order  rate  constants  (k0bs)  in  a  manner  that  systematically  accounts  for  the  absorbance  of 
colloidal  MnCL,  we  fit  all  of  our  kinetic  data  to  a  two-term  model: 


425  =[C;.e^]x^7  +  [C, 


,525 
MnO 9 


[6] 


1  Waldemer,  R.  H.;  Tratnyek,  P.G.  (2006)  Kinetics  of  contaminant  degradation  by  permanganate. 
Environ.  Sci.  Technol.  40(3):  1055-1061. 

2  Waldemer,  R.  H.;  Tratnyek,  P.G.  (2004)  The  efficient  determination  of  rate  constants  for 
oxidations  by  permanganate.  Proceedings  of  the  Fourth  International  Conference  on 
Remediation  of  Chlorinated  and  Recalcitrant  Compounds,  24-27  May  2004,  Monterey,  CA, 
Paper  2A-09. 

Waldemer,  R.  H.  (2004)  Determination  of  the  Rate  of  Contaminant  Degradation  by 
Permanganate:  Implications  for  In  Situ  Chemical  Oxidation  (ISCO).  M.S.  thesis,  OGI  School 
of  Science  and  Engineering,  Portland,  OR. 
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where  A  525  is  the  absorbance  at  525  nm,  C,  is  the  initial  amount  of  MnOf,  is  the 
absorptivity  of  MnOf,  and  £l^0i  is  the  absorptivity  of  Mn02.  We  treated  C,  and  as 
constants  and  simultaneously  fit  k0bs  and  by  least  squares  regression.  This  treatment  allows 
for  variability  in  ■>  which  is  expected  because  the  optical  properties  of  the  particles  will 
depend  on  their  size  distribution,  and  which  will  be  affected  by  the  rate  of  MnOf  reduction, 
which  will  depend  on  the  COC  that  is  being  oxidized.  The  fitted  values  we  obtained  for 
£mIo2  have  a  roughly  normal  distribution  (statistical  analysis  given  in  (Waldemer,  2004)). 
Overall,  equation  6  fit  the  experimental  data  very  well,  as  illustrated  in  Figure  2 A  for  a  fast 
reaction  (3-chlorophenol,  with  a  half-life  of  seconds)  and  Figure  2B  for  a  slow  reaction  (1,4- 
dioxane,  with  a  half-life  of  days).  Further  details  on  the  derivation  and  validation  of  equation  6 
can  be  found  in  (Waldemer,  2004;  Waldemer  and  Tratnyek,  2004). 


Figure  22:  Data  for  the  reactions  of  MnOf  with  3-chlorophenol  and  1 ,4-dioxane  (note  the 
difference  in  time  scales)  with  the  fit  to  a  simple  first-order  kinetic  model  (curve  in  red)  and 
Mn02-corrected  model  represented  by  equation  6  (curve  in  blue). 


The  k0bs  values  obtained  by  fitting  equation  6  to  experimental  data  were  used  to  determine 
second-order  rate  constants,  k".  For  COCs  with  low  solubility  (<  0.002  M)  or  reaction  rates  (half- 
lives  >1.5  days),  A:"  was  calculated  with  equation  7: 

*"=KJ\.coc]  [7] 

For  these  COCs,  the  treatment  effectively  assumes  the  order  of  reaction  with  respect  to  the 
COC  is  one.  The  COCs  for  which  k"  was  determined  this  way  were  MTBE,  toluene,  PCE,  picric 
acid,  and  2,4-dinitrophenol.  MTBE,  PCE,  and  toluene  have  been  shown  by  others  to  be  first- 
order  with  respect  to  the  COC  (Gardner,  1996;  Yan  and  Schwartz,  1999;  Huang  et  al.,  2001; 
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Damm  et  al.,  2002).  Picric  acid  and  2,4-dinitrophenol  were  assumed  to  be  first-order  with  respect 
to  COC  based  on  their  structural  similarity  to  the  mononitrophenols. 

For  COCs  with  sufficiently  high  solubilities  that  experiments  could  be  set  up  over  a 
significant  range  of  initial  COC  concentrations  ([COC]o),  we  plotted  k0bs  vs.  [COC]o,  determined 
that  the  plots  were  linear  (and  therefore  approximately  first  order  in  [COC]),  then  fit  the  data  by 
linear  regression,  and  assigned  the  fitted  slope  to  k".  Figures  3A  and  3B  illustrate  this  analysis 
for  3-chlorophenol  and  1,4-dioxane,  respectively. 


Figure  23:  Data  for  the  reactions  of  Mn04"  with  3-chlorophenol  and  1 ,4-dioxane:  /cobs  vs. 
concentration  of  COC  plots;  the  slope  of  the  regression  line  is  used  to  determine  k". 

In  addition  to  absorption  by  Mn02,  another  factor  that  can  complicate  the  determination  of 
k"  by  measuring  the  disappearance  of  MnOf  rather  than  the  COC  is  that  MnOf  can  be  consumed 
by  processes  other  than  the  primary  reaction  of  interest:  for  example,  the  autodecomposition  of 
MnOf  (even  in  pure  water  (Stewart,  1965))  and  reaction  with  daughter  products  derived  from 
the  initial  oxidation  of  the  COC.  Under  the  conditions  used  in  this  study,  however,  we  found  that 
autodecomposition  of  MnOf  was  negligible.  We  also  concluded  that  reaction  of  MnOf  with  any 
daughter  products  of  the  COC  oxidation  was  not  likely  to  be  significant  because  the  COC 
concentration  was  high — at  least  5-fold  (and  often  >  10-fold)  greater  than  the  MnOf — and 
therefore  would  be  the  dominant  reactant. 

Further  evidence  that  these  (or  other)  factors  did  not  significantly  bias  our  measurement  of 
k"  is  provided  by  the  strong  correlation  between  values  of  k”  obtained  from  this  study  and 
previously  published  values  of  k"  that  were  obtained  by  analyzing  decreasing  concentrations  of 
COC  in  the  presence  of  excess  MnOf  (Figure  4).  In  Figure  4,  perfect  agreement  between  our 
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data  and  literature  data  would  fall  on  a  line  with  slope  =  1  and  intercept  =  0,  which  is  shown  for 
comparison. 


-5  -4  -3  -2-10  1  2 

log  k"  with  excess  KMn04  (literature  data) 


Figure  24:  Correlation  between  /("obtained  from  pseudo-first-order  conditions  by  either 
analyzing  decreasing  concentrations  of  Mn04'  in  the  presence  of  an  excess  amount  of  COC 
(this  study)  or  analyzing  decreasing  concentrations  of  COC  in  the  presence  of  an  excess 
amount  of  Mn04"  (k"  data  obtained  from  the  literature).  Perfect  agreement  between  the  two  sets 
of  data  would  follow  the  dashed  line.  For  the  chlorinated  ethenes,  the  literature  values  of  k" 
used  in  this  Figure  are  the  values  Huang  et  al.  measured  at  25  °C  (Huang  et  al.,  2001)  because 
their  experimental  conditions  most  closely  matched  ours.  Multiple  literature  values  of  k"  for 
toluene  are  shown  because  the  only  literature  value  obtained  at  25  °C  was  also  obtained  with 
excess  COC  conditions.  For  MTBE  and  2,4-dichloro-phenol,  no  literature  data  from  25  °C  were 
available  for  this  comparison,  so  data  for  23  and  16  °C,  respectively,  were  used. 

Figure  5  summarizes  all  values  of  k"  that  were  measured  in  this  study  along  with  all  values 
of  k”  that  we  found  in  the  literature  for  Mn04‘  oxidation  of  potential  COCs  under  relevant 
conditions.  Inspection  of  the  data  in  Figure  5  shows  that  variability  in  k”  for  any  particular  COC 
is  generally  less  than  10-fold,  while  variability  among  closely  related  COCs  is  2-3  orders  of 
magnitude,  and  variability  among  families  of  unrelated  COCs  is  almost  7  orders  of  magnitude. 
The  variability  in  k”  for  any  particular  COC  is  satisfactory  considering  that  the  data  often  are 
from  different  laboratories  obtained  by  different  methods  under  slightly  different  conditions. 


12 


With  respect  to  the  variability  in  k"  within  and  across  families  of  COCs,  a  general 
observation  that  can  be  made  from  Figure  5  is  that  a  number  of  COCs  react  with  MnOf  nearly  as 
fast  or  faster  than  the  chlorinated  ethenes.  This  comparison  is  potentially  significant  because  the 
chlorinated  ethenes  have  been  remediated  successfully  with  MnOf  in  field-scale  applications  of 
ISCO  (Schnarr  et  al.,  1998;  Siegrist  et  al.,  2001),  and,  therefore,  it  follows  that  the  other  fast¬ 
reacting  COCs  shown  in  Figure  5  (e.g.,  TNT,  aldicarb,  dichlorvos,  many  phenols,  and  some 
PAHs)  are  possible  candidates  for  ISCO  with  MnOf.  Of  course,  having  established  that  k"  is 
suitably  large  for  these  COCs,  other  factors — such  as  the  effect  of  adsorption  to  soil  or  the 
potential  to  form  toxic  daughter  products — need  to  be  evaluated  to  establish  treatability.  Another 
useful  generalization  that  can  be  made  from  Figure  5  concerns  COCs  that  do  not  react  rapidly 
with  MnO/f.  The  families  of  COCs  with  values  of  k"  that  are  significantly  below  those  of  the 
chlorinated  ethenes  are  the  chlorinated  methanes  and  ethanes,  simple  aromatics  like  benzene  and 
toluene,  most  fuel  oxygenates  and  solvent  ethers,  and  some  explosives  (most  notably  RDX). 

As  part  of  this  study,  we  have  filled  many  of  the  priority  data  gaps  that  existed  previously, 
including  those  for  TNT,  1,4-dioxane,  some  pesticides,  and  many  phenols.  However,  because  our 
method  requires  that  COCs  have  at  least  moderate  solubility  (~0.5  mM),  we  were  not  able  to  fill 
some  data  gaps  for  significant,  low- solubility  compounds — most  notably  pentachlorophenol, 
dibenzodioxins,  polychlorinated  biphenyls,  and  some  of  the  more  insoluble  pesticides,  such  as 
l,l,l-trichloro-2,2-bis(p-chlorophenyl)ethane  (DDT),  and  dieldrin.  Most  of  these  COCs  are 
included  in  Figure  5 — even  though  no  k"  data  for  them  are  available  yet — so  that  the  Figure 
provides  a  balanced  perspective  on  current  coverage  of  the  whole  range  of  potential  COCs. 
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Figure  25:  Summary  of  second-order  rate  constants  (I k ")  for  Mn04'  and  compounds  in  all 
classes.  The  size  of  the  data  points  represent  the  temperature  at  which  the  data  were  collected; 
most  of  the  temperatures  for  the  literature  values  vary  between  16  °C  to  35  °C,  the  exception 
being  that  the  k"  reported  for  benzene  was  determined  at  70  °C.  The  pH  for  the  literature  values 
of  k"  ranges  from  4.6  to  8.0.  Compounds  labeled  in  red  text  are  too  insoluble  for  the  method  we 
used.  Experimental  conditions  for  the  data  points  from  this  study:  25  °C,  pH  7,  and  phosphate 
buffer  concentrations  between  50-  100  mM. 
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7.  Determining  K"  for  Activated  Hydrogen  Peroxide4 

An  exhaustive  literature  review  was  performed  to  summarize  existing  kinetic  data  on  the 
oxidation  of  NAPL  constituents  by  activated  hydrogen  peroxide  and  hydroxyl  radical.  As 
expected,  we  found  a  large  body  of  literature  on  the  oxidation  of  contaminants  by  classical 
Fenton’s  reagent.  However,  in  most  cases,  only  first  order  kinetics  are  reported.  In  this 
literature,  phenols  are  best  represented  with  24  papers  (e.g.,  Chen  and  Pignatello,  1997;  Esplugas 
et  ah,  2002)  and  oxygenates  are  least  represented  with  only  one  paper  (Burbano  et  ah,  2002). 

There  is  also  a  large  body  of  literature  reporting  k"  for  contaminant  degradation  by  *OH. 
However,  in  most  cases  the  *OH  is  generated  using  non-Fenton’s  methods,  such  as  gamma-  or 
pulse-radiolysis  (Farhataziz  and  Ross,  1977;  Buxton  et  ah,  1988),  photo-Fenton’s  (Haag  and 
Yao,  1992),  and  ozone  plus  hydrogen  peroxide  (Spanggord  et  ah,  2000).  In  this  literature,  all 
chemical  classes  are  well  represented  except  pesticides.  We  identified  only  about  ten  papers  that 
report  k”  of  contaminant  degradation  by  *OH  generated  by  classic  Fenton’s  reagent  (Haag  and 
Yao,  1992;  Tang  and  Huang,  1996b;  Tang  and  Huang,  1996a;  Gallard  et  ah,  1998;  Benitez  et  ah, 
2000;  Lindsey  and  Tarr,  2000;  Watts  et  ah,  2000;  De  Heredia  et  ah,  2001;  Gallard  and  De  Laat, 
2001).  These  ten  papers  include  representatives  of  all  chemical  classes  except  explosives. 

In  our  work  on  activated  hydrogen  peroxide,  we  set  out  to  design  experiments  that  would 
provide  (first-order)  kinetic  data  on  key  chlorinated  solvents  (priority  contaminants  for  which  k" 
are  lacking)  and  a  selective  trap  for  hydroxyl  radical  (benzoic  acid).  We  planned  to  use  the  latter 
data  to  normalize  the  former  data,  thereby  extracting  second-order  rate  constants  (as  illustrated  in 
Figure  1). 

Standard  experimental  conditions  were  choose  were  as  follows:  [H2O2]  =  8.8  mM,  [Fe2+] 

2_|_ 

=  0.45  mM  (as  FeSOzi),  and  [Fe  ]:[  H2O2]  molar  ratio  =  1:19.  All  experiments  used  house- 
distilled  water.  Reaction  mixtures  were  unbuffered  and  were  adjusted  to  an  initial  pH  of  3.0 
using  sulphuric  acid.  All  Fenton  experiments  are  being  performed  in  a  zero-headspace  reactor, 
adapted  from  (Chen  et  ah,  2001a;  Chen  et  ah,  2001b),  which  is  shown  in  Figure  6. 


4  We  have  no  publications  planned  from  this  part  of  the  project. 
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Figure  26.  Zero-headspace  reactor  used  for  to  measure  degradation  kinetics  by  iron-activated 
hydrogen  peroxide.  Sampling  is  performed  via  the  syringe. 


Experiments  were  performed  with  a  variety  of  chlorinated  solvents — including 
trichloroethene  (TCE)  and  1,1-dichloroethene  (DCE) — over  a  range  of  conditions.  As  an 
example,  the  conditions  used  for  TCE  are  presented  in  Table  1. 


Table  1.  Initial  conditions  in  the  Fenton’s  System  with  TCE. 


Exp 

[TCE] 

(mg/1) 

[H2o2] 

(mg/1) 

[Fe2+] 

(mg/1) 

[TCE] 

(mmol/1) 

[H202] 

(mmol/1) 

[Fe2+] 

(mmol/1) 

TCE/H202/Fe2+ 
molar  ratio 

1-2 

60 

300 

25 

0.46 

8.8 

0.45 

1:19:1 

3-4 

35 

300 

25 

0.29 

8.8 

0.45 

0.64:19:1 

5-6 

10 

300 

25 

0.08 

8.8 

0.45 

0.002:19:1 

The  results  we  obtained  with  each  chlorinated  solvent  are  illustrated  in  Figure  7  for  TCE 
and  DCE.  Clearly,  these  data  cannot  be  described  by  simple  pseudo  first-order  disappearance 
kinetics.  Using  a  more  complex  kinetic  model  that  we  developed,  we  were  able  to  fit  most 
aspects  of  these  data  (Figure  8). 
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Figure  27.  Disappearance  kinetics  for  chlorinated  solvents  in  a  homogenous  Fenton  system: 
(A)  TCE,  and  (B)  1,1  -DCE.  Solid  lines  illustrate  fits  to  pseudo  first-order  kinetics. 


Time  (min)  Time  (min) 

Figure  28.  Modeling  of  kinetics  for  chlorinated  solvents  in  a  homogenous  Fenton  system:  (A) 
TCE,  and  (B)  1,1-DCE.  Solid  lines  are  globally  fit  (across  all  sets  of  data  in  each  plot)  to  a 
subset  of  the  data. 

Similarly  complex  kinetics  were  obtained  with  our  hydroxyl  radical  trap  benzoic  acid  (not 
shown).  In  an  attempt  to  obtain  a  comprehensive  interpretation  and  model  for  these  data  we 
collaborated  with  colleagues  at  the  University  of  New  South  Wales  who  were  just  completing  a 
multiyear  effort  to  develop  a  fully-mechanistic  model  for  the  kinetics  of  Fenton  system  reactions 
(Duesterberg  and  Waite,  2006;  Duesterberg  and  Waite,  2007).  In  principle,  such  a  model  could 
be  used  to  extract  high  quality  second-order  rate  constants  for  the  reactions  of  *OH,  test  the 
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potential  significance  of  other  reactive  oxygen  species  (e.g.,  superoxide),  etc.  Unfortunately,  we 
were  not  able  to  get  this  model  to  fit  all  the  key  features  of  our  data.  We  expect  that  there  is  some 
inconsistency  in  the  data  and/or  the  model,  but  it  was  beyond  the  resources  of  this  project  to 
resolve  this.  Instead,  we  turned  out  focus  to  activated  persulfate,  for  which  kinetic  data  for 
environmentally-relevant  conditions  are  far  more  scarce. 

8.  Determining  k"  for  Activated  Persulfate’5 

Persulfate  (S2O8  ')  was  a  more  challenging  oxidant  to  work  with  than  permanganate  and 
similar  to  H2O2  in  that  most  of  the  observed  oxidation  of  contaminants  occurs  via  activation  to 
•OH  in  the  case  of  H2O2  and  SOzf”  in  the  case  of  S2O82".  The  two  general  ways  of  activating 
S2O8  '  are:  homolysis  of  the  peroxide  bond  using  heat  or  light  (equation  8)  and  an  oxidation- 
reduction  process  (analogous  to  the  Fenton  reaction)  with  electron  donors,  including  e'  (aq)  from 
radiolysis  of  water  or  low-valent  metals  (Mn+)  such  as  Fe2+  and  Ag+  (equation  9). 


S2082-  heatorlight  >2SO;~  [8] 

S2082'+Mn+ - >Mn+1+S0/-+S042-  [9] 

Once  sulfate  radical  is  generated,  the  following  reactions  have  been  speculated  to  occur: 

S0;-+H20<==>:H0'+H++S02-  [10] 

2  so;- — >s2o2-  [H] 

so;-  +  ho' — >hso;  ^ 

2  HO  - >H202  [13] 

H202 - »H20  +  0.5  02  [14] 

H202  +  S202" - >2  H++2S02"  +  02 


Determining  k”  for  persulfate  is  complicated  because  it  is  the  active  species  formed 
(SOD — not  the  persulfate  itself — for  which  k"  must  be  obtained.  Often,  to  obtain  k", 
experiments  are  designed  so  that  the  reaction  of  the  oxidant  (Ox)  with  the  COC  is  pseudo  first- 
order.  The  second-order  rate  law  can  be  written  in  two  ways: 


5  Waldemer,  R.  H.;  Tratnyek,  P.G.  (2007)  Sulfate  radical  (SO4*  )  oxidation  kinetics  of  organic 
contaminants  in  aquatic  media.  Environ.  Sci.  Technol.,  in  prep. 
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d[COC]  _fcn[coc][Ox] 
di 


[16] 


®^  =  A:"[COC][Ox]  [17] 

di 


Therefore,  the  experiment  can  be  designed  so  either  (/),  the  [COC]  or  (//),  the  [Ox]  is  in 
significant  excess  (usually  lOx  the  concentration  of  the  other)  and  is  effectively  constant 
throughout  the  reaction;  in  other  words,  pseudo-first-order  conditions. 

However,  a  third  option — the  one  employed  here — is  to  use  a  competition  kinetics  method, 
in  which  k”  is  measured  by  d[probe]/dt,  where  the  probe  is  a  compound  such  as  p- 
nitrosodimethylaniline  (PND A)  that  is  easy  to  analyze  (PNDA  absorbs  strongly  at  440  nm)  and 
is  highly  reactive  with  the  oxidant  (for  example,  k”  for  PNDA  and  SO4'  is  2.3  x  109  M'V1).  The 
experimental  set  up  for  this  method  consists  of  several  batch  reactions  with  varying 
concentrations  of  COC  and  identical  concentrations  of  PNDA  and  SO4'  .  In  batch  reactions  with 
higher  COC  concentrations,  the  COC  is  better  able  to  compete  with  the  PNDA  for  available 
SO4*  (Fig  9A).  The  data  can  be  linearized  so  k"  can  be  determined  (Fig  9B). 


Figure  29:  Typical  competition  kinetic  plots.  (A)  Change  in  absorbance  at  440  nm  (indicating 
change  in  concentration  of  PNDA  over  90  min);  (B)  Linear  version  of  1st  six  data  points,  used  to 
calculate  k". 

The  main  advantage  of  the  competition  kinetics  method  is  that  it  allows  k"  for  S04*  to  be 
determined  without  knowing  the  concentration  of  SO4’  in  the  reaction,  which  is  difficult  to  do 
without  the  use  of  a  stop- flow  apparatus,  often  an  expensive  piece  of  equipment.  The  potential 
pitfalls  of  this  method  include  the  possibility  that  PNDA  could  react  with  the  COC  or  its 
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daughter  products  in  addition  to  the  probe.  In  either  case,  simple  competition  kinetics  will  not 
apply,  and  the  S-shaped  curve  typical  of  competition  kinetics  will  not  be  obtained  (Figure  10).  A 
second  problem  would  arise  if  the  daughter  products  react  with  S2O82"  to  produce  additional 
SO/-;  again,  the  S-shaped  curve  would  not  be  obtained.  In  this  case,  at  low  COC  concentrations, 
the  PNDA  is  depleted  more  than  it  is  when  no  COC  is  present  due  to  the  additional  SO/  in 
solution,  but  at  high  COC  concentrations,  the  COC  is  able  to  effectively  compete  with  PNDA 
despite  the  additional  SO/  in  solution. 
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Figure  30:  Data  for  perchloroethene  (PCE)  does  not  conform  to  the  S-shaped  curve  typical  of 
competition  kinetics. 

Although  the  competition  method  has  several  possible  pitfalls,  it  must  be  remembered  that 
issues  arise  with  the  more  “traditional”  methods  of  measuring  k"  as  well.  In  Section  6,  we 
discussed  possible  problems  with  measuring  d[Ox]/dt:  namely,  the  possibility  that  the  oxidant 
could  react  with  daughter  products  or  could  be  consumed  by  processes  other  than  the  primary 
reaction  of  interest  (in  the  case  of  permanganate,  this  would  be  autodecomposition,  while  in  the 
case  of  SO//  this  would  be  the  reactions  described  by  equations  10-12).  In  the  case  of 
persulfate,  it  is  also  possible  for  the  daughter  products  to  react  with  S20g2'  and  produce 
additional  SO/  ,  which  would  make  the  value  of  k"  artificially  low.  Similarly,  when  obtaining 
k "  by  measuring  d[COC]/dt,  it  is  possible  for  the  COC  to  react  with  species  other  than  SO/  that 
are  generated  from  the  activation  of  S2O82"  (see  equations  10-15).  Furthermore,  for  a  study  such 
as  this  one,  where  values  of  k"  for  many  COCs  with  very  different  chemical  compositions  are 
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desired,  this  approach  is  not  practical,  because  a  different  analytical  approach  would  have  to  be 
developed  for  each  type  of  compound. 

In  light  of  all  the  possible  issues  with  all  methods  of  determining  k",  it  is  important  to 
validate  the  data  by  making  comparisons  among  multiple  methods  (as  we  did  for  permanganate 
in  Figure  4).  The  data  from  this  study  using  competition  kinetics  is  compared  with  data  from  the 
literature,  which  was  obtained  by  measuring  either  d[SCV  ]/dt  or  d[COC]/dt  (Figure  1 1).  The 
fact  that  the  data  (with  the  exception  of  MTBE)  land  on  a  one-to-one  line  is  indicative  that  all 
methods,  despite  their  potential  pitfalls,  are  appropriate  for  measuring  k".  Figure  1 1  also  shows 
that  in  the  experiments  outlined  for  this  study,  the  reactive  species  is  S(V~,  not  *OH. 
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Figure  31 :  The  values  of  k"  obtained  by  this  study  with  competition  kinetics  compared  with 
values  of  /c"for  ‘OH  and  S04’“  obtained  with  either  COC  or  S04‘“  in  excess  (literature  values). 


The  compilation  of  the  data  generated  by  competition  experiments  from  this  study  and  the 
rate  constants  obtained  by  the  literature  are  shown  in  Figures  12  (aromatic  compounds)  and  13 
(nonaromatic  compounds).  We  have  filled  key  data  gaps  for  toluene,  ethylbenzene,  2,4- 
dinitrotoluene,  2,6-dinitrotoluene,  1,1,1-trichloroethane,  and  1,1,2-trichloroethane.  Furthermore, 
we  have  added  data  for  MTBE,  t-butyl  alcohol,  chloroform,  and  ethanol.  The  chlorinated 
ethenes  are  a  noticeable  gap  in  the  data;  unfortunately,  when  the  competition  method  was  used 
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with  these  compounds,  S-shaped  curves  were  not  obtained,  indicating  that  k"  could  not  be  found 
using  this  probe. 

When  looking  at  Figures  12-13,  one  can  see  that  for  compounds  with  replicate  data 
generated  by  multiple  groups  the  agreement  is  almost  always  within  half  an  order  of  magnitude. 
This  is  especially  notable  when  considering  this  holds  true  even  when  SO4*  was  generated  in 
different  ways  (i.e.  narrowband  light,  broadband  light,  heat,  or  iron).  Thus,  as  expected,  it 
appears  that  the  method  of  persulfate  activation  does  not  greatly  affect  the  value  of  k”  for 
oxidation  of  COCs  by  SO4’  . 
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Figure  32.  Summary  of  second-order  rate  constants  (k")  for  S04'"  and  aromatic  compounds. 
The  big  circles  and  error  bars  represent  data  obtained  from  this  study,  small  circles  represent 
literature  values.  The  data  are  color  coded  by  the  method  used  to  activate  the  persulfate. 
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Figure  33.  Summary  of  second-order  rate  constants  (k")  for  S04‘~  and  aliphatic  compounds. 
The  big  circles  and  error  bars  represent  data  obtained  from  this  study,  small  circles  represent 
literature  values.  The  data  are  color  coded  by  the  method  used  to  activate  the  persulfate. 
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9.  Correlation  Analysis6 

One  of  the  goals  of  this  study  was  to  investigate  whether  quantitative  structure-activity 
relationships  (QSARs)  could  be  developed  for  oxidations  of  COCs  by  MnOf,  SO/-,  and  *OH. 
The  primary  motivation  for  developing  QSARs  was  to  provide  a  means  of  predicting  the  rate  of 
oxidation  for  compounds  that  have  not  been  studied.  A  further  benefit  of  developing  QSARs  is 
that  evidence  for  a  correlation  (or  lack  of  one)  between  chemical  structure  and  reactivity  with  an 
oxidant  can  provide  insight  on  the  mechanism  of  oxidation  with  COCs.  Because  outliers  often 
are  easily  identified  when  data  are  shown  in  this  format,  yet  another  benefit  of  QSARs  is  its 
usefulness  in  determining  consistency  among  the  data  (Tratnyek,  1998). 

In  order  to  develop  a  successful  QSAR  among  a  group  of  compounds,  the  compounds 
must  have  similar  reaction  mechanisms  or  structurally  similar  reaction  centers.  In  addition,  a 
descriptor  variable  that  relates  to  the  reaction  mechanism  or  reaction  center  must  be  identified. 
Descriptors  that  relate  to  the  energetics  of  reaction  include  the  half-wave  potential  (Ei/2), 
ionization  potential  (IP),  and  the  energy  of  the  highest  occupied  molecular  orbital  (EHomo)- 
These  may  be  preferred  descriptors  for  oxidation  reactions  (Tratnyek,  1998;  Canonica  and 
Tratnyek  Paul,  2003). 

The  descriptor  that  gave  the  best  correlation  for  oxidation  of  COCs  with  SO/-  was  Ehomo 
(Figure  14),  but  the  correlation  is  to  complex  to  fit  with  a  simple  univariate  QSAR  model. 
Regretfully,  we  did  not  have  the  resources  to  explore  more  complex  (nonlinear)  models  that 
might  explain  the  data,  we  do  not  yet  have  a  QSAR  to  report  that  can  be  used  to  predict  values  of 
k"  for  the  oxidation  reactions  by  SO/-  for  COCs  that  have  not  been  previously  studied. 

However,  the  fact  that  the  data  set  used  for  the  correlation  includes  all  the  compounds  shown  in 
Figures  12  and  13,  and  the  basic  trend — that  compounds  with  low  values  of  EHomo  are  least 
reactive  with  SO/- — is  upheld  for  the  entire  set,  lends  confidence  that  for  most  COCs,  those 
with  high  values  of  EHomo  will  be  most  reactive  with  SO/-.  Also  shown  in  Figure  14  is  that 
aromatic  compounds  tend  to  have  the  highest  reactivity  with  SOzf”,  while  nonaromatic 
compounds  without  double  or  triple  bonds  tend  to  have  lower  reactivity  with  SO/-.  This 
observation  gives  insight  into  the  mechanism  of  these  oxidation  reactions,  implying  that  in 
addition  to  energetics,  the  electrophilic  nature  of  S(V_  also  plays  a  role  in  the  oxidation  of 
COCs. 


6  Waldemer,  Rachel.  H.;  Kim,  Jae-Hyoun;  Tratynek,  Paul  G.  (2007)  Correlation  analysis  of 
kinetic  data  for  aqueous  oxidants.  Journal  of  Chemical  Information  and  Modeling',  in  prep. 


25 


10  — I 


9- 

8- 

7- 


O 

c/) 

o)  6  — 


5- 


4- 


3- 


i —  i 

-12.5  -12.0 


•  aromatic 

•  nonaromatics  w/at  least  one  multiple  bond 

•  nonaromatics  w/all  single  bonds 


-11.5  -11.0  -10.5  -10.0  -9.5  -9.0  -8.5 

Ehomo  (PM5/H20) 


Figure  34.  Correlation  of  /("for  oxidation  reactions  of  COCs  by  S04‘  to  EHomo- 

Unlike  S04*_,  oxidation  by  Mn04'  of  a  large  group  of  COCs  containing  different  chemical 
subgroups  is  not  described  by  a  significant  correlation  to  either  EHomo  or  IP.  Another  descriptor 
that  was  explored  is  Eq ap:  the  energy  difference  between  the  HOMO  and  the  lowest  unoccupied 
molecular  orbital  (LUMO).  Eqap  has  been  used  as  a  measure  of  the  relative  stability  of  a 
compound  toward  chemical  reaction — compounds  with  larger  gaps  are  thought  to  have  lower 
reactivity  (Karelson,  2000;  Miehr  et  al.,  2004).  However,  the  stability  of  the  COC  as  measured 
by  Eoap  also  does  not  appear  to  predict  the  reactivity  with  Mn04~.  Interestingly,  the  descriptor 
that  seems  to  correlate  best  with  reactivity  of  Mn04"  among  compounds  of  different  chemical 
classes  is  the  reactivity  of  these  compounds  with  ozone  (Figure  15).  It  is  therefore  likely  that 
there  are  other  effects  in  addition  to  electron  transfer  (perhaps  steric  effects)  that  are  involved  in 
the  reaction  of  these  compounds  with  Mn04\  If  the  unknown  effects  also  affect  the  reactivity  of 
compounds  with  ozone,  these  effects  would  be  incorporated  in  the  rate  of  oxidation  by  ozone. 
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Figure  35.  Scatter  plot  matrix  of  log  k"  of  all  compounds  obtained  from  experiments  or  literature. 
All  k"  values  were  obtained  at  a  pH  between  4.0  and  8.0,  at  temperatures  between  20°C  and 
30°C.  Ehomo  and  EGapwere  calculated  with  CAChe  molecular  modeling  software.  IP  values 
were  obtained  from  (Watanabe  et  al.;  Kobayashi  and  Nagakura,  1974;  Kobayashi  and 
Nagakura,  1975).  Ozone  rate  constants  were  obtained  from  (Hoigne  and  Bader,  1983b;  Hoigne 
and  Bader,  1983a;  Hoigne  et  al.,  1985). 


Although  a  satisfactory  correlation  could  not  be  found  to  describe  oxidation  by  MnO/f  for 
the  entire  dataset  of  COCs,  a  meaningful  correlation  was  found  for  oxidation  of  the  chlorinated 
ethenes  (Figure  16).  Because  chlorine  is  an  electron-withdrawing  substituent,  and  the  reactivity 
of  the  chlorinated  ethenes  decreases  with  increasing  number  of  chlorines,  MnOzf  appears  to  be 
behaving  as  an  electrophile  in  the  initial  rate-determining  step,  and  it  is  likely  that  the  electron 
density  in  the  pi  bonds  of  the  chlorinated  ethylenes,  and  not  energetics,  is  controlling  the 
reaction. 
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Figure  36.  Matrix  of  scatter  plots  of  log  k" for  Mn04‘  vs.  the  number  of  chlorine  atoms  attached 
to  the  double  bond  and  vs.  log  k'Tor  03.  (Ozone  data  from  (Hoigne  and  Bader,  1983a)). 
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In  addition  to  the  correlation  of  the  rate  constants  for  MnO/f  and  the  chlorinated  ethenes 
to  the  number  of  chlorine  atoms  on  the  double  bond,  there  is  also  a  satisfactory  correlation  with 
k"  for  the  oxidation  of  these  compounds  by  ozone.  Direct  ozonation  of  unsaturated  compounds 
in  low-pH  aqueous  media  (the  media  used  in  the  referenced  ozone  experiments)  is  generally 
agreed  to  occur  through  the  formation  of  ozonide,  another  five-membered  cyclic  intermediate 
(Larson  and  Weber,  1994);  therefore,  it  likely  incorporates  many  of  the  rate-determining  effects 
that  are  not  described  by  the  other  descriptors.  It  is  not  surprising,  then,  that  the  rates  of 
oxidation  between  ozone  and  MnOT  correlate  so  well — not  only  is  there  a  high  value  of  R2  for 
the  correlation  (0.86),  but  it  is  clear  that  the  correlation  would  be  even  better  by  removing  the 
cis-DCE  outlier.  That  the  reaction  rate  of  oxidation  of  cis-DCE  by  MnOT  is  over-predicted  by 
the  rate  of  oxidation  by  ozone  is  likely  because  ozonide  is  less  bulky  than  the  cyclic 
hypomanganate  diester  formed  by  MnOT  oxidation,  and  is  therefore  not  as  greatly  affected  by 
steric  hindrance  from  the  close  proximity  of  the  chlorine  atoms  on  cis-DCE. 

Given  that  reactivity  with  ozone  was  the  best  descriptor  for  the  entire  dataset  of  kinetic 
data  with  MnOT,  upon  gathering  data  for  S(V_  and  "OH  we  included  it  in  the  cross-correlation 
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matrix.  As  can  be  seen  in  Figure  15,  the  only  correlations  with  values  of  R  higher  than  0.5  are 
between  SO/-  and  "OH,  SO/-  and  O3,  "OH  and  O3,  and  MnOzf  and  O3.  That  SOzf-  and  "OH 
correlate  with  each  other  but  not  MnOT  is  likely  indicative  of  the  non-specific  nature  of  the  first 
two  oxidants  as  compared  to  MnOf.  That  all  three  of  the  aforementioned  oxidants  correlate  to 
ozone  may  then  seem  surprising;  however,  this  observation  is  easily  explained  by  the  fact  that 
ozone  can  oxidize  compounds  directly,  in  the  manner  similar  to  MnOf  that  was  described 
earlier,  or  indirectly,  through  production  of  "OH  (Larson  and  Weber,  1994).  It  is  likely  that  the 
compounds  that  correlate  well  with  SO4*  and  'OH  are  oxidized  by  the  indirect  route,  and  the 
compounds  that  correlate  well  with  M11O4  are  oxidized  directly.  The  oxidation  of  chloroform 
(CHCI3)  illustrates  this  point  well.  Chloroform  is  clearly  an  outlier  in  the  correlations  of 'OH 
and  S04"  with  ozone,  but  not  in  the  correlation  of  Mn04_,  which  may  indicate  that  chloroform  is 
oxidized  by  ozone  directly. 

10.  IscoKin  database7 

In  addition  to  our  publications,  we  created  a  database  containing  our  collection  of  k"  data 
that  is  available  on  the  Internet  at  http://cgr.ebs.ogi.edu/iscokin.  Access  is  unrestricted  and  we 
plan  to  keep  the  database  online  for  the  foreseeable  future.  A  screen  shot  of  the  search  interface 
is  shown  in  Figure  17. 

Unlike  many  online  databases,  which  do  not  provide  for  automatic  updating  as  new  data 
become  available,  our  database  allows  registered  users  to  input  new  values  of  k".  This  feature 
makes  it  possible  for  the  ISCOKIN  database  to  remain  a  current  and  complete  resource  for 
practitioners  of  ISCO  and  researchers  interested  in  correlation  analysis. 


n 

http://cgr.ebs.ogi.edu/iscokin 
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A  database  of  second  order  rate  constants  for  contaminant  remediation  using  chemical  oxidants. 
For  background  and  other  details,  see  the  IscoKin  database  Help  file. 


Figure  37.  Screen  shot  of  the  initial  search  window  for  the  IscoKin  database. 
(http://cgr.ebs.ogi.edu/iscokin/search.php) 

11.  Oxidation  of  Chlorinated  Ethenes  by  Heat-Activated  Persulfate8 

A  potential  advantage  of  S2O8  is  its  stability  before  activation,  which  may  allow 
delivery  of  the  oxidant  to  contamination  in  hard  to  reach  places.  To  take  full  advantage  of  this 
property,  the  best  method  of  activating  S2O82  may  be  the  remote,  localized,  and  directed  heating 
provided  by  some  of  the  technologies  used  in  ISTR.  (For  a  review  of  ISTR  heating  technologies, 
see  (Fountain,  1998).)  This  approach  to  activating  S2O8  contrasts  favorably  with  activation 
methods  that  involve  mixing  reagents  before  injection  (e.g.,  S2O8  with  chelated  iron)  because 
the  latter  inevitably  result  in  reaction  of  some  of  the  oxidant  before  it  reaches  the  contaminated 
zone.  In  the  case  of  S2O82  activation  with  chelated  iron,  we  expect  that  some  of  the  SO4’ 
produced  will  react  with  excess  Fe2+  ( k "  =  3-9.9  x  108  NT's"1  (Heckel  et  al.,  1966;  McElroy  and 
Waygood,  1990))  thereby  lowering  the  concentration  of  SO/-  that  is  available  to  degrade 
contaminants.  This  effect  may  account  for  the  observation  that  iron  activation  of  S2O82  is  not 
always  effective  at  degrading  some  contaminants  that  are  degraded  with  S2O82  activated  by 
other  methods  (Block  et  ah,  2004). 

As  mentioned  earlier,  once  S2082-  is  activated  (equations  8  or  9),  the  resulting  S(V_ 
initiates  a  chain  of  reactions  involving  other  radicals  and  oxidants  (equations  10-15),  some  of 
which  are  potentially  reactive  intermediates,  such  as  "OH  and  peroxymonosulfate  (HSO5  ).  The 


8  Waldemer,  R.  H.;  Tratnyek,  P.G.;  Johnson,  Richard  L.;  Nurmi,  James  T.  (2007)  Oxidation  of 
chlorinated  ethenes  by  heat  activated  persulfate:  Kinetics  and  products.  Environ.  Sci.  Technol,. 
31(3),  1010-1015. 
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variety  of  intermediate  oxidants  generated  in  activated  S2O8  systems  complicates  the  kinetics 
of  contaminant  oxidation.  In  principle,  the  contaminant  disappearance  should  be  describable  with 
a  pseudo-first-order  rate  constant  that  is  the  sum  of  second-order  terms  for  each  oxidant: 

=  Ko,-  [«V  ]  +  k"„o[°H]  +  KA-  [>A‘  ]  +  C„,  [other]  [1 8] 

where  k"  represents  the  second-order  rate  constants  for  the  reaction  of  the  contaminant  with  each 
reactive  intermediate.  Under  most  conditions,  the  dominant  term  in  equation  18  is  presumed  to 
be  the  one  involving  SO4*  (Kolthoff  et  al.,  1951;  Peyton,  1993;  Couttenye  et  al.,  2002),  which  is 
why  knowing  k”  for  the  oxidation  of  SO4’  and  COCs  is  valuable.  However,  it  is  not  known  how 
the  relative  significance  of  these  terms  varies  with  system  parameters,  such  as  temperature.  Of 
course,  k"  for  all  the  terms  in  equation  18  will  increase  with  temperature,  resulting  in  faster 
oxidation  of  contaminants  if  the  oxidant  concentrations  are  not  decreased.  The  latter,  however, 
can  not  be  assumed  because  the  concentrations  of  S(V_  and  other  oxidants  shown  in  equation  18 
are  the  net  result  of  many  reactions  (equations  8-15,  the  reactions  of  each  oxidant  with  the 
contaminants  and  their  byproducts,  etc.),  each  of  which  will  have  its  own  dependence  of  rate  on 
temperature.  An  additional  uncertainty  is  the  effect  of  temperature  on  the  relative  rates  of 
contaminant  degradation  by  chemical  oxidants  other  than  S2O8  and  background  degradation 
processes  such  as  hydrolysis.  This  study  addresses  these  issues  for  the  chlorinated  ethenes  by 
determining  the  temperature  dependence  of  the  kinetics  and  products  of  their  oxidation  by 
S2O8  ,  and  comparing  these  results  to  calculations  made  from  previously  published  data  on  the 
degradation  of  the  chlorinated  ethenes  by  permanganate  and  hydrolysis. 

As  shown  in  Figure  18  for  PCE,  oxidation  of  the  chlorinated  ethenes  with  thermally- 
activated  Na2S208  was  pseudo-first-order  for  1-2  half-lives  at  every  temperature  studied.  In  some 
cases,  data  were  collected  for  more  than  two  half-lives  and  the  later  time-points  exhibited  the 
distinctive  downward  curvature  shown  in  Figure  19.  The  fact  that  this  curvature  can  also  be  seen 
in  previously  published  data  for  TCE  and  1,1,1  -trichloroethane  oxidation  by  S2O8  (Liang  et  al., 
2003)  suggests  that  the  effect  is  not  peculiar  to  the  conditions  of  this  study.  The  cause  of  this 
curvature  was  not  explored  in  this  or  the  previous  studies,  but  it  could  be  due  to  the  combined 
reactivity  of  SO4’  and  other  reactive  species  that  form  as  the  reaction  progresses,  such  as  the 
organic  radicals  and  chlorine  radical  (Cl")  that  may  be  by-products  of  oxidation  of  chlorinated 
ethenes.  To  avoid  having  to  model  this  phenomenon,  we  fit  only  the  initial  linear  portion  of  the 
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data  (1-2  half-lives)  using  the  pseudo-first-order  model.  The  resulting  values  of  k0bs  are  given  in 
Table  2  for  the  four  chlorinated  ethenes  and  five  temperatures  studied. 


Figure  38.  Pseudo-first-order  disappearance  of  PCE  at  five  temperatures.  Experimental 
conditions:  unbuffered,  4.5  *  10"5  M  PCE,  and  4.5  *  10'4  M  Na2S208.  Triplicate  experiments  at 
each  temperature.  The  control  was  done  without  persulfate. 
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Figure  39.  Concentration  over  extended  time  periods,  showing  non-pseudo  first-order  behavior 
at  longer  experimental  times.  Experimental  conditions:  Unbuffered  Dl  water,  70  °C,  10:1  molar 
ratio  of  Na2S208  and  chlorinated  ethenes,  concentration  of  PCE  =  0.0014  M,  concentration  of 
the  other  chlorinated  ethenes  =  0.002  M. 
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Table  2.  kobs  and  corresponding  activation  parameters  for  all  4  chlorinated  ethenes 


Compound 

30°C 

40  °C 

50  °C 

60  °C 

70  °C 

In  A 

(s'1) 

Ea 

(kJmol'1) 

0.023  hr'1 

0.95 

0.24  hr'1 

0.99 

0.0087  min'1 

0.94 

0.041  min'1 

0.98 

0.11  min'1 

0.97 

0.045  mM 

PCE 

0.027  hr'1 

0.97 

0.25  hr'1 

0.96 

0.0070  min'1 

0.91 

0.030  min"1 

0.98 

0.080  min'1 

0.99 

0.031  hr1 

0.97 

0.24  hr'1 

0.96 

0.0070  min'1 

0.97 

0.034  min'1 

0.94 

0.083  min"1 

0.97 

0.060  hr1 

0.99 

— 

— 

— 

— 

101  ±4 

0.055  hr1 

0.99 

— 

— 

— 

— 

"  29  ±  2 

0.0014  M 

PCE 

0.038  hr1 

0.99 

— 

— 

— 

0.148  min'1 

0.99 

0.040  hr'1 

0.99 

— 

— 

— 

0.145  min'1 

0.98 

— 

— 

— 

— 

0.1 15  min"1 

0.98 

0.002  M 

TCE 

0.027  hr1 

0.98 

0.073  hr'1 

0.94 

0.0046  min'1 

0.96 

0.021 1  min"1 

0.99 

0.061  min'1 

0.99 

0.031  hr'1 

0.97 

0.087  hr'1 

0.97 

0.0045  min'1 

0.97 

0.0184  min'1 

0.99 

0.067  min'1 

0.99 

31  ±  1 

108  ±3 

0.026  hr'1 

0.99 

0.076  hr'1 

0.93 

0.0051  min"1 

0.96 

0.0179  min'1 

0.99 

0.069  min"1 

0.9 

0.002  M 

cis-DCE 

0.0033  hr1 

0.93 

0.023  hr'1 

0.96 

0.0017  mm1 

0.94 

0.008  min'1 

0.91 

0.056  min"1 

0.99 

0.0034  hr1 

0.96 

0.022  hr'1 

0.99 

0.0016  min'1 

0.98 

0.0054  min'1 

0.99 

0.068  min"1 

0.98 

43  ±2 

144  ±5 

0.0032  hr1 

0.95 

0.026  hr'1 

0.99 

— 

0.005  min"1 

0.80 

0.057  min"1 

0.99 

0.002  M 

trans-DCE 

0.0048  hr1 

0.97 

0.033  hr'1 

0.99 

0.0030  min"1 

0.98 

0.018  min'1 

0.98 

0.050  min"1 

0.99 

0.0060  hr1 

0.96 

0.038  hr'1 

0.99 

0.0030  min"1 

0.99 

0.017  min'1 

0.98 

0.058  min"1 

0.99 

43  ±  1 

141  ±2 

0.0056  hr1 

0.97 

0.032  hr'1 

0.97 

0.017  min'1 

0.97 

0.059  min'1 

0.98 
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For  all  four  chlorinated  ethenes,  In  k0bs  decreased  linearly  with  1/T  (Figure  20)  and,  therefore, 
were  fit  to  an  Arrhenius-type  model: 

\nkobs=\nA-EA/RT  [19] 

where  A  is  the  pre-exponential  factor,  EA  is  the  apparent  activation  energy,  R  is  the  universal  gas 
constant,  and  T  is  the  absolute  temperature.  The  fitting  parameters  from  these  calculations  are 
presented  in  Table  3.  Both  In  A  and  EA  follow  the  pattern  PCE  <  TCE  <  trans-DCE  <  cis-DCE;  in 
fact,  In  A  and  EA  correlate  well  with  R  =  0.99  (Waldemer  et  al.,  2007),  a  phenomenon  that  is 
fairly  common  but  of  uncertain  significance  (Philibert,  2006).  A  consequence  of  this  is  that  kQbs 
also  follows  the  pattern  PCE  >  TCE  >  trans-DCE  >  cis-DCE  and  correlates  with  In  A  and  EA 
(Waldemer  et  al.,  2007).  Additionally,  In  A,  EA  ,  and  k0bs  all  correlate  fairly  well  with  the  energy 
of  the  highest  occupied  molecular  orbital  ( EHomo )  for  the  chlorinated  ethenes.  This  is  to  be 
expected  because  oxidation — regardless  of  the  exact  mechanism — generally  involves  a  shift  of 
electron  density  from  the  HOMO  of  the  reductant  (in  this  case,  the  chlorinated  ethenes)  to  an 
electron  acceptor  (in  this  case,  SO/f  )(Canonica  and  Tratnyek,  2003). 
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Figure  40  Arrhenius  plots  for  the  chlorinated  ethenes.  Experimental  conditions:  unbuffered, 
10:1  molar  ratio  of  Na2S208  to  chlorinated  ethenes.  Initial  concentrations  are  given  in  the 
legend. 


Table  3:  Arrhenius  parameters  for  chlorinated  ethenes. 


Compound 

In  A  (s'1)* 

Ea  (kJ/mol)* 

N 

R2 

PCE 

29  ±2 

101  ±4 

25 

0.96 

TCE 

31  ±  1 

108  ±3 

15 

0.99 

cis-DCE 

44  ±2 

144  ±5 

14 

0.98 

trans-DCE 

43  ±  1 

141  ±2 

14 

0.99 

Uncertainties  are  one  standard  deviation  derived  from  fitting  equation  19  to  the  data  in  Figure 
20. 


Our  value  of  EA  for  TCE  is  within  10%  of  the  value  reported  previously  (97.7  kJ  mol'1) 
by  Liang  et  al.  (Liang  et  al.,  2003).  The  only  other  data  on  EA  for  oxidation  of  chlorinated 
ethenes  by  heat-activated  S208  are  from  Huang  et  al.  (Huang  et  al.,  2005).  These  values  (PCE: 
46.4,  TCE:  60.7,  cis-DCE:  41.4,  trans-DCE:  49.8,  1,1-dichloroethene:  64.9,  and  vinyl  chloride: 
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56.9  kJ  mol'1)  are  significantly  lower  (2-3  fold)  than  the  values  reported  here  (and  by  Liang  et  al. 
(Liang  et  al.,  2003)),  and  the  two  sets  of  data  do  not  correlate  well  (not  shown).  A  likely 
explanation  for  this  disagreement  lies  in  the  fact  that  Huang  et  al.  measured  degradation  kinetics 
in  a  mixture  of  59  contaminants  (Huang  et  al.,  2005).  In  such  a  situation,  all  of  the 
contaminants — and  possibly  their  degradation  intermediates  and  products — can  compete  for 
SOf  (and  any  other  reactive  species).  If  the  various  reactions  have  different  activation  energies, 
some  contaminants  will  have  proportionally  higher  reaction  rates  at  the  same  temperature  than 
others,  and  the  contaminants  with  higher  activation  energies  will  become  more  effective 
competitors  at  higher  temperatures.  This  competition  effect  could  explain  why  the  apparent 
activation  energies  for  the  chlorinated  ethenes  are  lower  in  pseudo-first-order  reactions  with 
many  other  species  present  than  when  they  are  the  only  reductant  in  solution. 

Little  has  been  reported  on  the  products  of  oxidation  of  the  chlorinated  ethenes  by  S20g2  , 
and  there  is  no  information  on  how  the  reaction  pathways  and  products  are  affected  by 
temperature.  Anticipating  some  degree  of  mineralization,  we  measured  CL  present  at  the  end  of 
experiments  performed  with  the  four  chlorinated  ethenes  at  70  °C.  For  this  purpose,  we  chose  to 
run  each  experiment  for  a  period  anticipated  to  encompass  4  half-lives  of  contaminant 
degradation  (-94%),  which  was  estimated  from  the  average  of  the  pseudo-first-order  rate 
constants  listed  in  Table  2.  Overall,  CL  recoveries  were  consistent  with  complete  dechlorination 
of  80-90%  of  each  chlorinated  ethene  (data  not  shown).  The  only  comparable  literature  data  are 
for  TCE,  and  they  include  CL  recoveries  of  75%  at  40°C,  80%  at  50°C,  and  85%  at  60°C  (Liang 
et  al.,  2003). 

For  systems  containing  SOf  ,  however,  the  interpretation  of  Cl  data  is  complicated  by 
the  possibility  that  they  react  according  to  equation  20  (k”  =  2.47  x  10  M's'  (Huie  and  Clifton, 
1990)).  Then,  the  resulting  chlorine  radical,  Cf,  can  react  with  additional  CL,  establishing  an 
equilibrium  with  Cf'  according  to  equation  21  (k"=  8  x  109  M'V1  (Huie  and  Clifton,  1990),  K 
=  1.4  x  105  M'1  (Buxton  et  al.,  2000)). 


so;~  +  cr — 

->cr+so42- 

[20] 

cr+ci 

cr2~ 

[21] 

Equation  21  suggests  that  the  predominant  chlorine  radical  species  at  the  initiation  of  the  Cl 
recovery  experiments  was  Cf ,  when  the  initial  Cl  concentration  was  very  low  (below  detection), 
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and  that  CF*  should  have  become  predominant  by  the  end  of  each  experiment  (because 
degradation  of  the  chlorinated  ethenes  produced  final  CF  concentrations  2-5  mM,  which  should 
shift  equation  21  to  favor  Cl2*  ).  One  consequence  of  equations  20-21  could  be  depletion  of  Cl  , 
but  of  greater  interest  is  the  possibility  that  Cl*  and  CF*  react  with  the  chlorinated  ethenes  or 
their  intermediate  degradation  products  to  produce  more  highly  chlorinated  products,  as  has  been 
reported  previously  for  chlorinated  phenols  (Anipsitakis  et  al.,  2006).  The  kinetic  data  necessary 
to  evaluate  how  favorable  chlorination  of  ethenes  might  be  do  not  appear  to  be  available,  but  we 
did  detect  trace  quantities  of  hexachloroethane  as  intermediates  in  the  oxidation  of  TCE  and 
PCE. 

The  only  chlorinated  compounds  that  were  identified  as  major  intermediates  were  cis- 
DCE  during  trans-DCE  oxidation  and  trans-DCE  during  cis-DCE  oxidation.  This  evidence  for 
racemization  is  shown  in  Figure  21  for  70  °C,  but  the  effect  was  seen  at  all  temperatures. 
Racemization  presumably  occurs  after  formation  of  a  single -bonded  intermediate,  which  might 
be  the  radical  cation  formed  by  electron-transfer  to  S(V_  (Minisci  et  al.,  1983),  or  the  S(V_- 
radical  adduct,  formed  by  asymmetric  addition  of  S(V_  to  the  double  bond.  Evidence  for  the 
latter  is  strong,  with  many  SO4*  -radical  adducts  having  been  identified  during  reactions  of  SO4* 
with  alkenes  by  electron  spin  resonance  (Chawla  and  Fessenden,  1975;  Koltzenburg  et  al.,  1982; 
Davies  and  Gilbert,  1984).  It  is  less  clear,  however,  what  factors  control  the  elimination  of  S(V_ 
to  form  the  racemized  DCE’s. 


Figure  41 .  Concentration  of  cis-DCE  and  trans-DCE  vs.  time  for  A)  the  reaction  of  cis-DCE  with 
persulfate  and  B)  the  reaction  of  trans-DCE  with  persulfate.  Experimental  conditions: 
unbuffered,  triplicate  experiments,  0.002  M  dichloroethenes,  0.02  M  Na2S208,  70  °C. 
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To  evaluate  the  potential  benefits  of  coupling  ISCO  and  ISTR  for  chlorinated  ethenes,  it 
is  necessary  to  consider  how  temperature  affects  a  variety  of  contaminant  degradation  processes. 
For  hydrolysis  and  oxidation  by  permanganate,  it  is  straightforward  to  calculate  the  respective 
values  of  k0bs  vs.  temperature  (up  to  100  °C)  using  equation  19  previously  published  Arrhenius 
parameters  (Jeffers  et  al.,  1989;  Huang  et  al.,  2001),  and  appropriate  assumptions  regarding  the 
dose  of  permanganate  (100-40,000  mg  L'1,  0.0007  M  to  0.28  M  (Siegrist  et  al.,  2001))  and 
groundwater  pH.  For  PCE,  we  show  the  results  of  these  calculations  in  Figure  22,  with  the 
thickness  of  the  bands  reflecting  the  relevant  ranges  of  permanganate  concentration  and  pH.  For 
the  other  chlorinated  ethenes,  we  used  only  a  representative  permanganate  concentration  and  pH, 
and  the  results  are  shown  in  Figure  23.  Over  the  whole  range  of  conditions  considered,  oxidation 
of  the  chlorinated  ethenes  is  much  faster  than  hydrolysis,  although  the  opposite  would  have  been 
the  case  for  some  of  the  more  highly  chlorinated  ethanes.  Additional  factors  must  be  considered 
for  temperatures  over  100  °C  (Costanza  et  al.,  2005). 
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1.16  x  1014  yrs 

0  20  40  60  80  100 

Temp  (°C) 

Figure  42.  Comparison  of  three  degradation  processes  for  PCE  as  a  function  of  temperature. 
The  hydrolysis  band  was  calculated  using  Arrhenius  parameters  obtained  from  (Jeffers  et  al., 
1989).  The  bottom  of  the  pink  band  corresponds  to  100  mg/L  (0.63  mM)  permanganate  and  the 
to  of  this  band  corresponds  to  40,000  mg/L  (253  mM),  calculated  using  Arrhenius  parameters 
obtained  from  (Huang  et  al.,  2001).  The  light  blue  band  corresponds  to  Arrhenius  parameters 
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obtained  from  experiments  with  0.45  mM  S208/0.045  mM  PCE.  The  blue  dots  correspond  to 
experimental  data  for  varying  concentrations  of  S208  while  keeping  the  PCE  concentration  at 
0.045  mM.  The  three  concentrations  represented  by  the  three  increasingly  dark  blue  circles  are: 
4.5  mM  S208,  45  mM  S208,  and  450  mM  S208.  The  dashed  blue  line  allows  the  reader  to  follow 
the  likely  Arrhenius  progression  of  the  450  mM  data  at  all  temperatures. 


Figure  43.  Comparison  of  three  degradation  processes  as  a  function  of  temperature  for  TCE, 
cis-DCE,  and  trans-DCE.  The  hydrolysis  band  was  calculated  using  Arrhenius  parameters 
obtained  from  (Jeffers  et  al.,  1989).  The  permanganate  band  was  calculated  using  Arrhenius 
parameters  obtained  from  (Huang  et  al.,  2001).  The  persulfate  band  was  calculated  using 
Arrhenius  parameters  obtained  from  experiments  with  20  mM  Na2S208  and  2  mM  of  the 
respective  chlorinated  ethenes 


For  contaminant  oxidation  by  activated  S20g2  ,  it  is  more  difficult  to  develop  a  general 
description  of  k0bs  vs.  T  because  there  are  many  reactions  that  can  contribute  to  the  concentration 
of  the  reactive  intermediate  SO4*  (equations  8-15,  etc.),  each  with  its  own  dependence  on 
temperature.  (Note  that  calculating  the  temperature  dependence  of  contaminant  oxidation  by  the 
Fenton  reaction  would  pose  a  similar  challenge.)  To  overcome  this,  we  used  the  EA  for  PCE  that 
was  determined  in  this  study  (Table  3)  to  fix  the  slope  of  the  temperature  dependence  and 

2 _ 

estimated  the  intercept  (In  A)  from  replicate  experiments  done  at  initial  concentrations  of  S2Og 
ranging  from  107  -  107,000  mg/L  (0.00045  M  to  0.45  M).  This  amounts  to  assuming  that  the 
apparent  Ea  is  unaffected  by  S20g2  concentration.  The  result,  shown  in  Figure  22,  reveals  that 
the  temperature  at  which  the  rate  of  PCE  oxidation  by  heat-activated  S2Og  becomes  greater 
than  oxidation  by  permanganate  is  approximately  40  °C;  and,  by  100  °C,  heat-activated  S20g2- 
oxidizes  PCE  about  400  times  faster  than  permanganate.  In  contrast,  the  other  chlorinated 
ethenes  (Figure  23)  are  so  much  more  reactive  with  permanganate  that  oxidation  with  activated 
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S2O82  is  slower  even  at  100  °C.  The  relative  reactivity  of  PCE  with  these  two  oxidants  is 

consistent  with  the  relative  reactivity  of  the  chlorinated  ethenes  with  each  oxidant:  PCE  >  TCE  > 

2 _ 

cis-DCE  >  trans-DCE  with  activated  SAT  (Figure  20)  and  the  reverse  with  permanganate 
(Waldemer  and  Tratnyek,  2006). 

While  the  comparisons  facilitated  by  Figures  22  and  23  should  be  qualitatively  accurate, 
the  absolute  rates  we  calculated  may  differ  from  those  observed  in  the  field  for  a  variety  of 
reasons.  Prominent  among  the  factors  that  suggest  slower  degradation  rates  under  in  situ 
conditions  is  the  inefficiency  of  mixing  of  aqueous-phase  oxidants  with  aqueous-  and 
nonaqueous-phase  contaminants  because  of  aquifer  heterogeneity  (Seol  et  al.,  2003),  limited 
dispersion  (especially  vertical  mixing  (Johnson  and  Pankow,  1992))  and  displacement  of  non¬ 
sorbing  contaminants  by  the  injected  fluids.  This  limitation  might  be  ameliorated  by  coupling 
ISCO  with  ISTR,  because  ISTR  should  encourage  mixing  of  the  oxidant  plume  with  the 
contamination  (by  increasing  flow  within  the  heated  zone  due  both  to  reduced  water  viscosity 
and  increased  buoyancy  of  the  heated  water).  Another  reason  to  expect  slower  oxidation  rates  in 
the  field  is  that  “natural  oxidant  demand”  will  suppress  SO/-  concentrations  (Brown  and 
Robinson,  2004)  leaving  less  SO4'  to  oxidize  contaminants.  Here,  again,  coupling  ISCO  with 
ISTR  may  be  advantageous,  because  S2O8  will  be  relatively  stable  in  unheated  zones  and  heat- 
activation  can  be  focused  on  contaminated  zones.  Additional  factors  remain  to  be  investigated, 
such  as  thermal  effects  on  geochemical  processes  that  might,  in  turn,  influence  contaminant 
oxidation  rates. 

12.  Effect  of  groundwater  constituents  on  ISCO  kinetics9,10 

In  general,  we  expect  that  k"  as  defined  throughout  this  project  will  depend  on  pH, 
temperature,  and  perhaps  ionic  strength,  but  not  the  concentrations  of  other  groundwater 
constituents;  whereas  k0bs  will  also  depend  on  the  concentration  of  the  oxidant  and  therefore  be 


9  Waldemer,  R.  H.;  Tratnyek,  P.G.;  Johnson,  Richard  L.;  Nurmi,  James  T.  (2007)  Oxidation  of 
chlorinated  ethenes  by  heat  activated  persulfate:  Kinetics  and  products.  Environ.  Sci.  Technol,. 
31(3),  1010-1015. 

10  Waldemer,  R.  H.  (2004)  Determination  of  the  Rate  of  Contaminant  Degradation  by 
Permanganate:  Implications  for  In  Situ  Chemical  Oxidation  (ISCO).  M.S.  thesis,  OGI  School 
of  Science  and  Engineering,  Portland,  OR. 
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impacted  by  any  processes  that  compete  with  the  COC  for  the  oxidant  (Hoigne,  1990;  Tratnyek 
and  Macalady,  2000).  The  latter  effect  should  be  negligible  for  MnOf  because  it  reacts  directly 
with  COCs,  and  this  is  consistent  with  the  observation  that  the  rate  of  TCE  oxidation  by  MnOf  is 
the  same  in  “hard”  tap  water  and  in  deionized  water  (Vella  and  Veronda,  1993).  In  contrast, 
activated  persulfate  and  activated  H2O2  react  with  COCs  via  reactive  radicals  that  can  be 
scavenged  by  dissolved  carbonate  and  other  species,  which  in  turn  makes  kQbs  for  these  treatment 
processes  sensitive  to  the  concentrations  of  such  groundwater  constituents  (2005). 

To  verify  that  bicarbonate  and  other  common  groundwater  constituents  do  not  have  a 
significant  effect  on  rates  of  COC  oxidation  by  MnOf,  we  measured  k0bs  for  oxidation  of  TCE 
over  a  range  of  concentrations  for  the  major  groundwater  constituents.  Bicarbonate 
concentrations  in  the  range  of  50  to  400  mg/L  often  occur  in  groundwater,  with  concentrations 
above  1000  mg/L  possible;  nitrate  concentrations  in  groundwater  are  mostly  below  20  mg/L,  but 
groundwater  contaminated  with  fertilizers  can  have  concentrations  above  1000  mg/L;  and  sulfate 
concentrations  in  groundwater  tend  to  range  from  under  15  mg/L  to  about  500  mg/L  (Matthess, 
1982).  Therefore,  we  chose  to  study  the  effect  of  0  to  1000  mg/L  nitrate  (0-15.9  mM),  sulfate 
(0-10.4  mM),  and  bicarbonate  (0-16.4  mM)  on  the  oxidation  of  TCE  (1.0  mM)  by  MnOf  (0.1 
mM).  As  expected,  the  rates  ( k0bs  and  k”)  were  not  affected  by  any  of  the  above  groundwater 
constituents  over  the  range  of  concentration  studied  (data  shown  in  (Waldemer,  2004)). 

As  mentioned  earlier,  unlike  oxidation  reactions  with  MnOf,  we  did  expect 
environmental  concentrations  of  carbonate  species  to  affect  kabs  for  oxidation  by  activated 
S2O82  .  To  determine  the  degree  of  this  effect  on  heat- activated  S20g2-,  PCE  was  used  as  the 
model  compound  along  with  100,  300,  and  500  mg  L'1  of  bicarbonate  at  30,  50,  and  70  °C.  The 
data  were  fit  to  pseudo  first-order  kinetics — as  shown  in  figure  22  for  30  °C — and  the  resulting 
values  of  kabs  are  given  in  Table  4.  At  all  temperatures,  bicarbonate  inhibited  the  oxidation  of 
PCE  by  activated  S2O82  (figure  23).  The  inhibitory  effect  of  bicarbonate  appears  to  be  steeper 
than  a  simple  exponential  relationship  and  equivalent  at  all  three  temperatures,  but  we  did  not 
attempt  to  model  this  behavior. 
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Figure  44.  Effect  of  bicarbonate  on  the  oxidation  of  PCE  by  persulfate.  Experimental  conditions: 
30  °C,  4.5  x  10'5  M  PCE,  4.5  x  10'4  M  Na2S208,  two  replicates.  The  control  did  not  contain 
persulfate  or  bicarbonate. 


The  inhibitory  effect  of  bicarbonate  probably  reflects  competition  between  the  PCE  and 
bicarbonate  for  SCV”,  but  changes  in  pH  could  also  be  contributing.  For  all  the  experiments 
shown  in  Figure  23,  the  final  pH  decreased — possibly  due  to  the  release  of  protons  when  S(V_ 
reacts  with  water  (equation  10) — and  the  degree  of  pH  decrease  was  less  pronounced  with 
increased  bicarbonate  concentration  (Figure  22).  However,  even  with  bicarbonate  present,  the 
oxidation  of  PCE  still  follows  the  Arrhenius  model  (shown  in  (Waldemer  et  al.,  2007))  and  the 
activation  energies  are  not  statistically  different  with  bicarbonate  or  without  it  (Table  3). 
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Table  4:  k0 bs  and  corresponding  activation  parameters  for  PCE  with  different 
concentrations  of  bicarbonate 


Bicarbonate 

(mg  L'1) 

30°C 

50°C 

70°C 

In  A 

(s'1) 

Ea 

(kJ  mol'1) 

0 

See  Table  1 

See  Table  1 

See  Table  1 

29  ±2 

101  ±4 

100 

0.018 

0.30 

0.84 

25  ±2 

94  ±5 

300 

0.013 

0.12 

1.62 

26  ±5 

97  ±  12 

500 

0.011 

0.15 

1.32 

26  ±  1 

96  ±2 

Figure  45.  Effect  of  bicarbonate  concentration  on  kobs  for  PCE  oxidation  by  persulfate  activated 
at  three  temperatures.  4.5  x  10'5  M  PCE,  4.5  *  10'4  M  Na2S208. 
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